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I. INTRODUCTION 



Stationary vortex loops are often discussed in the literature in vario us con texts 



Torn the models of condensed matter physics 111 | 




P, J4l| J124 Jl37j, [l46|, [l47|, 



U, [id], [177 



ra ngin g 



149| . to high energy physics and cosmology 29|, 45|, 46|, (48|, 16^ . 1771 ] . Such objects 
could be quite interesting physically and might be responsible for a plenty of important 
physical phenomena, starting from the structure of quantum superfiuids 5J] and Bose- 
Einstein condensed alkali g ases |ll7l ] to the baryon asymmetry of the Universe, dark matter, 
and the galaxy formation 1651 ] . These intriguing physical aspects of vortex loops occupy 



therefore most of the discussions in the hterature, while much less attention is usually 
payed to their mathematical existence. In most considerations the existence of such loops is 
discussed only qualitatively, using plausibility arguments, and the problem of constructing 
the corresponding field theory solutions is rarely addressed, such that almost no solutions 
of this type are explicitly known. 

The physical arguments usually invoked to justify the loop existence arise within the 
effective macroscopic description of vortices [l|, [l29|. Viewed from large distances, their 
internal structure can be neglected and the vortices can be effectively desc ribed as elastic 
thin ropes or, if they carry currents, as thin wires [32], 33|, 3^, |35|], 122| . This suggests 
the following engineering procedure: to 'cut out' a finite piece of the vortex, then 'twist' 
it several times and finally 'glue' its ends. The resulting loop will be stabilized against 
contraction by the potential energy of the twist deformation, that is by the positive pressure 
contributed by the current and stresses |43|], 42|, 86|, 130||. Equally, instead of twisting the 
rope, one can first 'make' a loop and then 'spin it up' thus giving to it an angular momentum 
to stabilize it against shrinking [47|, 48|. Depending on whether they are spinning or not, 
such loops are often called in the literature vortons and knots (also springs), respectively. 

These macroscopic arguments are suggestive. However, they cannot guarantee the exis- 
tence of loops as stationary field theory objects, since they do not take into account all the 
field degrees of freedom that could be essential. Speaking more rigorously, one can prepare 
twisted or spinning loop configurations as the initial data for the fields. However, nothing 
guarantees that these data will evolve to non-trivial equilibrium field theory objects, since, 
up to few exceptions, the typical field theory models under consideration do not have a 
topological energy bound for static loop solitons. For spinning loops the situation can be 
better, since the angular momentum can provide an additional stabilization of the system. 
However, nothing excludes the possibility of radiative energy leakages from spinning loops. 
Since there is typically a current circulating along them, which means an accelerated motion 
of charges, it seems plausible that spinning loops should radiate, in which case they would 
not be stationary but at best only quasistationary. In fact, loop formation is generically 
observed in dynamical vortex network simulations (although for currentless vortices), and 
these loops indeed radiate rapidly away all their energy (see e.g. 26|). 

The lack of explicit solutions renders the situation even more controversial. In fact, 
although there are explicit examples of loop solitons in global field theory, which means 



containing only scalar fields with global internal symmetries, almost nothing is known about 
such solutions in gauge field theory. Vortons, for example, were initially proposed more than 
20 years a go in the context of the local U(l) xU(l) theory of superconducting cosmic strings 



of Witten (l77j ]. However, not a single field theory solution of this type has been obtained 
up to now. Search for static knot solitons in the Ginzburg-Landau type gauge field theory 
models has already 30 years of history, the result being always negative. All this casts some 
doubts on the existence of stationary, non-radiating vortex loops in physically interesting 
gauge field theory models. However, rigorous no-go arguments are not known either. As a 
result, there are no solid arguments neither for nor against the existence of such solutions. 
Interestingly, exact solutions describing ring-type objects are known in curved space. 
These are black rings in the multidimensional generalizations of General Relativity - spinning 



toroidal black holes (see [58j for a review). It seems therefore that Einstein's field equations, 
notoriously known for their complexity, are easier to solve than the non-linear equations 
describing loops made of interacting gauge and scalar fields in Minkowski space. 

It should be emphasized that the problem here is not related to the dynamical stability 
of these loops. In order to be stable or unstable they should first of all exist as stationary 
field theory solutions. The problem is related to the very existence of such solutions. One 
should be able to decide whether they exist or not, which is a matter of principle, and 
this is the main issue that we address in this paper. Therefore, when talking about loops 
stabilized by some forces, we shall mean the force balance that makes possible their existence, 
and not the stability of these loops with respect to all dynamical perturbations. If they 
exist, their dynamical stability should be analyzed separately, but we shall call solitons all 
localized, globally regular, finite energy field theory solutions, irrespectively of their stability 
properties. It should also be stressed that we insist on the stationarity condition for the 
solitons, which implies the absence of radiation. Quasistationary loops which radiate slowly 
and live long enough could also be physically interesting, but we are only interested in loops 
which could live infinitely longtime, at least in classical theory. 

Trying to clarify the situation, we review in what follows the known field theory solu- 
tions in Minkowski space in 3 + 1 dimensions describing stationary ring objects. We shall 
divide them in two groups, depending on whether they do have or do not have an angular 
momentum. Solutions without angular momentum are typically static, that is they do not 
depend on time and also do not typically have an electric field. We start by describing 
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the static knot solitons stabilized by the Hopf charge in a global field theory model. Then 
we review the attempts to generalize these solutions within gauge field theory, with the 
conclusion that some additional constraints on the gauge field are necessary, since fixing 
only the Hopf charge does not seem to be sufficient to stabilize the system in this case. We 
then consider two known examples of static ring solitons in gauge field theory: the anoma- 
lous solitons stabilized by the Chern-Simons number, and the non-Abelian monopole and 
sphaleron rings. 

Passing to the rotating solutions, we first of all discuss the issue of how the presence of 
an angular momentum, associated to some internal motions in the system, can be reconciled 
with the absence of radiation which is normally generated by these motions. One possibility 
for this is to consider time-independent solitons in theories with local internal symmetries. 
They will not radiate, but it can be shown for a number of important field theory models 
that the angular momentum vanishes in this case. Another possibility arises in systems with 
global symmetries, where one can consider non-manifestly stationaij and axisymmetric fields 
containing spinning phases. Such fields could have a non-zero angular momentum, but the 
absence of radiation is not automatically guaranteed in this case. The general conclusion 
is that the existence of non-radiating spinning solitons, although not impossible, seems to 
be rather restricted. Such solutions seem to exist only in some quite special field theory 
models, while generic spinning field systems should radiate. 

Nevertheless, non-radiating spinning solitons in Minkowski space in 3+1 dimensions exist, 
and we review below all known examples: these are the spinning Q-balls, their twisted and 
gauged generalizations, spinning Skyrmions and rotating monopole-antimonopole pairs. In 
addition, there are also vortons, and below we present for the first time numerical solutions 
of elliptic equations describing vortons in the global field theory limit. Although global 
vortons have been studied before by different methods, our construction shows that they 
indeed exist as stationary, non-radiating field theory objects. 

We also discuss stationary ring solitons in non-relativistic physics. Surprisingly, it turns 
out that the relativistic vortons can be mapped to the 'skyrmion' solutions of the Gross- 
Pitaevskii equation in the theory of Bose-Einstein condensation. In addition, it turns out 
that Q-balls can describe light pulses in media with non-linear refraction - 'light bullets'. 
Finally, for the sake of completeness, we discuss also the moving vortex rings stabilized by 
the Magnus force in continuous media theories, such as in the superfiuid helium and in 



ferromagnetics. 

In our numerical calculations we use an elliptic PDE solver with which we have managed 
to reproduce most of the solutions we describe, as well as obtain a number of new results 
presented below for the first time. The latter include the explicit vorton solutions, spinning 
twisted Q-balls, spinning gauged Q-balls, as well as the 'Saturn', 'hoop' and bi-ring solutions 
for the interacting Q-balls. We put the main emphasise on describing how the solutions are 
constructed and not to their physical applications, so that our approach is just the opposite 
and therefore complementary to the one generally adopted in the existing literature. As 
a result, we outline the current status of the ring soliton existence problem - within the 
numerical approach. Giving mathematically rigorous existence proofs is an issue that should 
be analyzed separately. It seems that our global vorton solutions could be generalized in 
the context of gauge field theory. The natural problem to attack would then be to obtain 
vortons in the electroweak sector of Standard Model. 

In this text the signature of the Minkowski spacetime metric g^^ is chosen to be 
(+, — , — , — ), the spacetime coordinates are denoted by x^ = {x^, x^) = {t, x) with k = 1,2,3. 
All physical quantities discussed below, including fields, coordinates, coupling constants and 
conserved quantities are dimensionless. 

II. KNOT SOLITONS 

Let us first consider solutions with zero angular momentum stabilized by their intrinsic 
deformations. We shall start by discussing the famous example of static knotted solitons in 
a non-linear sigma model. Since this is the best known and also in some sense canonical 
example of knot solitons, we shall describe it in some detail. We shall then review the status 
of gauge field theory generalizations of these solutions. 

A. Faddeev-Skyrme model 

More than 30 years ago Faddeev introduced a field theory consisting of a non-linear 0(3) 



sigma model augmented by adding a Skyrme-type term 6l|], 62|. This theory can also be 



obtained by a consistent truncation of the 0(4) Skyrme model (see Sec JVB l! below). Its 



dynamical variables are three scalar fields n = n"" = (n^,n^,n^) constraint by the condition 

3 

n ■ n = V^n^n" = 1, 

a=l 

so that they span a two-sphere S^. The Lagrangian density of the theory is 



1 



£[n] = ,^^^ {d,n ■ d^n - T.^^n (2-1) 



where 

1 1 

- eabcn^d^n^d^n" = - 

The Lagrangian field equations read 



^/.^ = - eabcn^'d^n^d^n^ = 77 n ■ (9^n x d^u). (2.2) 



d^d^'n + d^J^^^'in x d^n) = (n ■ a^S^n) n . (2.3) 

In the static limit the energy of the system is 

^[n] = ^f ( {dkuf + {J^.kf) rf^x =E2 + E,. (2.4) 

Under scale transformations, x -^ Ax, one has E2 — * AE2 and E4 —* E^^/ K. The energy will 
therefore be stationary for A = 1 if only the virial relation holds, 

E2 = E4. (2.5) 

This shows that the four derivative term E4 is necessary, since otherwise the virial relation 
would require that E2 = 0, thus rulin g o ut all non-trivial static solutions - in agreement 



with the Hobart-Derrick theorem 92], 52 1 



Any static field n(x) defines a map M'^ -^ S"^. Since for finite energy configurations n 
should approach a constant value for |x| -^ 00, all points at infinity of M^ map to one point 
on S*^. Using the global 0(3)-symmetry of the theory one can choose this point to be the 
north pole of the 5*^, 

lim n(x) = noo = (0,0,1). (2.6) 

x|— >oo 

Notice that this condition leaves a residual 0(2) symmetry of global rotations around the 
third axis in the internal space, 

n^ + in^ -^ {n^ + m\''', n^ -^ n\ (2.7) 

The position of the field configuration described by n(x) can now be defined as the set of 
points where the field is as far as possible from the vacuum value, that is the preimage of the 
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point — Hoo antipodal to the vacuum +noo- This preimage forms a closed loop (or collection 
of loops) called position curve. Solitons in the theory can therefore be viewed as string-like 
objects, stabilized by their topological charge to be defined below. 

At the intuitive level, these solitons can be viewed as closed loops made of twisted vortices. 
Specifically, the theory admits solutions describing straight vortices that can be parametrized 
in cylindrical coordinates {p, z , (pj as n^ = cos 0(p), n^+in'^ = sin Q{p) 6*^^+*"*^, where n & Z 



is the vortex winding number lll| . The phase pz + mp thus changes both along and around 
the vortex and the vector n rotates around the third internal direction as one moves along 
the vortex, which can be interpreted as twisting of the vortex. It seems plausible that a loop 
made of a piece of length L of such a twisted vortex, where pL = 27rm, could be stabilized 
by the potential energy of the twist deformation. For such a twisted loop the phase increases 
by 27m after a revolution around the vortex core and by 27rm as one travels along the loop, 
where m G Z is the number of twists. If the loop is homeomorphic to a circle, then the 
product nm gives the value of its topological invariant: the Hopf charge. 

B. Hopf charge 

The condition ( 12. 6p allows one to view the infinity of M.^ as one point, thus effectively 
replacing M^ by its one-point compactification S^. Any smooth field configuration can 
therefore be viewed as a map 

n(x) : S^ -^ S\ (2.8) 

Any such map can be characterized by the topological charge A^[n] G 7r3(S'^) = Z known as 
the Hopf invariant. This invariant has a simple interpretation. The preimage of a generic 
point on the target 5^ is a closed loop. If a field has Hopf number A^ then the two loops 
consisting of preimages of two generic distinct points on S"^ will be linked exactly A^ times 
(see FigUl). 

Although there is no local formula for A^[n] in terms of n, one can give a non-local 
expression as follows. The 2-form ^ = | ^ikdx'^ A dx^ defined by Eg. (12.21) is closed, dJF = 0, 
and since the second cohomology group of S^ is trivial, H{S^) = 0, there globally exists a 
vector potential A = Akdx^ such that JF = dA. The Hopf index can then be expressed as 

Ar[n] = ^ I ^ijkAJ'jk c^'x. (2.9) 



10 




FIG. 1: Preimages of any two points rii and n2 on the target space S'^ are two loops. The number 
of mutual linking of these two loops is the Hopf charge A^[n] of the map n(x). Here the case of 
A^ = 1 is schematically shown. 



For any smooth, finite energy field configuration n(x) this integral is integer- valued 119 |. 

It can be shown t hat t he maximal symmetry of n(x) compatible with a non-vanishing 
Hopf charge is 0(2) lll| . It follows that spherically symmetric fields are topologically 
trivial. However, axially symmetric fields can have any value of the Hopf charge. Using 
cylindrical coordinates such fields can be parametrized as llll ] 



n^ + iv? 



^i{mip—nip) 



sin 9, 



n 



cos 6 



(2.10) 



where n,??! G Z and Q,ip are functions of p,z. Since n 
vanish for r 



rioo asymptotically, 9 should 



^/^T 



oo. The regularity at the 2;-axis requires for m ^ that 



should vanish also there. As a result, one has n = rioo both at the 2;-axis and at infinity, 
that is at the contour C shown in Figj2l Next, one assumes that n = — noo on a circle 
S around the 2;-axis which is linked to C as shown in Figl2] (more generally, one can have 
n = —rioo on several circles around the z-axis). The phase function ip is supposed to increase 
by 27r after one revolution along C. Since cos 6 interpolates between —1 and 1 on every 
trajectory from S to C, it follows that surfaces of constant 6 are homeomorphic to tori. 

The preimage of the point —noo consists of m copies of the circle 5*. The preimage of noo 
consists of n copies of the contour C. These two preimages are therefore linked mn times. 

One can also compute JF = dA according to (12. 2p . from where one finds 



A = n cos^ — dtp + m sin^ — dip 



(2.11) 



so that Al\T = nm cos^ y sin Q di/j A dO A dip. Inserting this to (12. 9p gives the Hopf charge 



iV[n] 



mn. 



(2.12) 



Following Ref. 158| . we shall call the fields given by the ansatz (12.101) A. 
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FIG. 2: Knot topology: the complex phase of the field in Eq. (j2.10p winds along two orthogonal 
directions: along S and along the contour C consisting of the z-axis and a semi-circle whose radius 
expands to infinity. A similar winding of phases is found for other systems to be discussed below: 
for vortons, skyrmions, and twisted Q-balls. 

Let us consider two explicit examples of the A^n field. Let us introduce toroidal coordi- 
nates {u, V, if} such that 



p = (R/r) sinhu, z = (R/r) sinw, 



(2.13) 



where r = coshw — cost" with u G [0, oo), v = [0, 27r). The correct boundary conditions for 
the field will then be achieved by choosing in the ansatz (12.101) 



e = e{u), tp = v, 



(2.14) 



where 6(0) = 0, e(cx)) = vr. 

Another useful parametrization is achieved by expressing n°- in terms of its complex 
stereographic projection coordinate 



W 



n^ + irt^ 



(2.15) 



The values VT = 0, oo correspond to the vacuum and to the position curve of the soliton, 
respectively. Passing to spherical coordinates {r.'d.Lp} one introduces 



(p = cosx(^) + ismx{r) cos^, a = sinx(r) sin-i^e"^, 



(2.16) 
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such that |0p + IxP = 1, where x(0) = vr, x(oo) = 0. A particular case of the field fl2.10p is 
then obtained by setting 

C. Topological bound and the N = 1,2 hopfions 

The following inequality for the energy (12 ■4p and Hopf charge f l2.9p has been established 



by Vakulenko and Kapitanski 163 1 



E[n]>c\N[u\f\ (2.18) 



where c = (3/16)'^/^ 111| . Its derivation is non-trivial and proceeds via considering a 



sequence of Sobolev inequalities. It is worth noting that a fraction al po wer of the topological 



charge occurs in this topological bound, whose value is optimal 119l |. On the ot her h and 



171 that 



it seems that the value of c can be improved, that is increased. Ward conjectures 
the bound holds for c = 1, which has not been proven but is compatible with all the data 
available. The existence of this bound shows that smooth fields attaining it, if exist, describe 
topologically stable solitons. Constructing them implies minimizing the energy (12. 4p with 
fixed Hopf charge (12. 9p . Such minimum energy configurations are sometimes called in the 
literature Hopf solitons or hopfions, and we shall call them fundamental or ground state 
hopfions if they have the least possible energy for a given N . 

Hopfions have been constructed for the first time for the lowest two values of the Hopf 

n 

charge, A^ = 1,2 by Gladikowski and Hellmund [79| and almost simultaneously (although 
somewhat more qualitatively) by Faddeev and Niemi [6J]. Gladikowski and Hellmund used 
the axial ansatz (I2.10p expressed in toroidal coordinates, with 6 = 6(m, f) and ip = v + 
ipo{u, v), where 9(0, v) = 0, 9(oo, v) = n. Assuming the functions 6(u, v) and ^o(^) ^0 to be 
periodic in v, they discretized the variables u, v and numerically minimized the discretized 
expression for the energy with respect to the lattice cite values of Q,ipo- They found a 
smooth minimum energy configuration of the An type for A^ = 1, while for A^ = 2 they 
obtained two solutions, A21 and A12, the latter being more energetic than the former. 

We have verified the results of Gladikowski and Hellmund by integrating the field equa- 
tions (12.30 in the static, axially symmetric sector. Using the axial ansatz (I2.10p . the az- 
imuthal variable ip decouples, and the equations reduce to two coupled PDE's for Q{p,z), 
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ip{p,z). Unfortunately, these equations are rather comphcated and it is not possible to re- 
duce them to ODE's by further separating variables, via passing to toroidal coordinates, say. 
In fact, we are unaware of any attempts to solve these differential equations. We applied 
therefore our numerical method described below in Sec JVIl to integrate them, and we have 
succeeded in constructing the first two fundamental hopfions. An and A21. For the A^ = 1 
solution the energy density is maximal at the origin and the energy density isosurfaces are 
squashed spheres, while for the N = 2 solution they have toroidal structure (see FigjS]). For 





FIG. 3: The energy isosurfaces for the A^ = 1 (left) and N = 2 (right) fundamental hopfions. 



the solutions energies E^- we obtained the values Ei = 1.22, E2 = 2.00. These feat ures agree 



with the results of Gladikowski and Hellmund and with those of Refs. 



90|. In particular, the same values for the energy are quoted by Ward 



m^n 



m 



89|, 



Ward also proposes a sim ple analytic approximation of the solutions based on the rational 
ansatz formula (I2.17P 172| . For n = m = 1 this formula gives 

iy = -f^e- (2.19) 

where /(r) = rcotx(r). Here the constant phase factor has been introduced to account for 
the residual 0(2) symmetry (12.71) . Minimizing the energy with respect F{r) it turns out 
that choosing 

/(r) = 0.453 (r - 0.878){r^ + 0.705r + 1.415) (2.20) 



gives for the energy a value which is less than 1% above the true minimum, Ei = 1.22 172 |. 
Eqs. (l2.19l) . (l2.20p provide therefore a good analytic approximation for the A^ = 1 hopfion. 
They show, in particular, that the position curve of the hopfion is a circle of radius R = 0.848 
and that for large r the field shows a dipole type behaviour. 



Tt + m^ ^ 2W ^ a -^. 



f2.211 
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Ward also suggests the moduli space approximation in which the | A^| = 1 hopfion is viewed as 
an oriented circle (see FigH]). There are six continuous moduli parameters: three coordinates 
of the circle center, two angles determining the position of the circle axis, and also the overall 
phase. Choosing arbitrarily a direction along the axis (shown by the vertical arrow in FigJl]), 
there are two possible orientations corresponding to the sign of the Hopf charge, changing 
which is achieved hj W ^ W*. Although for one hopfion the phase is not important, the 
relative phases of several hopfions determine their interactions. 

Ward conjectures that well-separated hopfions interact as dipoles with the maximal at- 
traction/repulsion when they are parallel/antiparallel, respectively, since the like charges 
attract in a scalar field theory. He verifies this conjecture numerically, and then numerically 
relaxes a field configuration corresponding to two mutually attracting hopfions. He discovers 
two possible outcomes of this process. If the two hopfions are initially located in one plane 
then they approach each other till the two circles merge to one thereby forming the A2,i 
hopfion with N = 2. If they are initially oriented along the same line then they approach 
each other but do not merge even in the energy minimum, where they remain separated by 
a finite distance. This corresponds to the Ai^2 hopfion, which is more energetic but locally 
stable. 

The A2.1 hopfion can be approximated by 

(x + iy^e"^ 



W 



f = 0.23(r - 1.27)(r + 0.44)(r + 0.16)(H - 2.15r + 5.09) (2.22) 



/ + il.55 zr' 

whose energy is 1.5 % ab ove the true minimum, and there is also a similar approximation 
for the A12 hopfion 172 1. 

N=2 



N=l 



N=-l 






2,1 



1,2 



FIG. 4: Schematic representation of the hopfions as oriented circles 
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D. Unknots, links and knots 



Similarly to the A^ = 1, 2 solutions, hopfions can be constructed within the axial ansatz 
f l2.10p also for higher values of n, m. However, for N > 2 they will not generically correspond 
to the global energy minima. This can be understood if we remember that the Hopf charge 
measures the amount of the twist deformation. Twisting an elastic rod shows that if it is 
twisted too much then the loop made of it will not be planar, since it will find it energetically 
favorable to bend toward the third direction. It is therefore expected that the ground state 
hopfions for higher values of A^ will not be axially symmetric planar loops, but 3D loops, 
generically without any symmetries. 



The interest towards this issue was largely stirred by the work of Faddeev and Niemi 64 1 



651], who conjectured that higher A^ hopfions should be not just closed lines but knotted 
closed lines, with the degree of knotedness expressed in terms of the Hopf charge. In other 
words, they conjectured that there could be a field theory realization of s table knots - the 
idea that had been put forward by Lord Kelvin in the 19-th century 162l | but never found 
an actual realization. This knot conjecture of Faddeev and Niemi had a large resonance and 
several groups had started large scale numerical simulations to look for knotted hopfions. 
Such solutions have indeed been found, although not with quite the same properties as had 
been originally predicted in Ref. 6J], 65| . 

The first, .eally asto„ishi,.g set of results has been reported bv Battye and SuteUffe Q, 
[151 ] , who managed to construct hopfions up to A^ = 8 and found the first non-trivial knot - 
the trefoil knot - for N = 7. Similar analyses h ave b e en t he n in dependently carried out by 
Hietarinta and Salo 89| , 90| and also by Ward 171| , 1721 ] , 176| . All groups performed the 
full 3D energy minimization starting from an initial field configuration with a given A^. Var- 
ious initial configurations were used, as for example the ones given by the rational ansatz 
fl2.17p . supplemented by non-axially symmetric perturbations to break the exact toroidal 
symmetry. The value of A^ being constant during the relaxation, the numerical iterations 
were found to converge to non-trivial energy minima, whose structure was sometimes com- 
pletely different from that of the initial configuration (an online animation of the relaxation 
process is available in 59|]). Several local energy minima typically exist for a given A^, some- 
times with almost the same energy, so that it was not always easy to know whether the 
minima obtained were local or global. Different initial configurations were therefore tried 
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to see if the minimum energy configurations could be reproduced in a different way. As a 
result, it appears that the global energy minima have now been identified and cro ss-ch ecked 



up to A^ = 7 J90|], |l58l |. after which the analysis has been extended up to A^ = 16 [158 
properties of the solutions can be summarized as follows. 



The 



For A^ = 1,2 these are the toroidal hopfions of Gladikowski and Hellmund 79|, An 
and A21. Although initially obtained within the constrained, axially symmetric relaxation 



the fu 



1 3D energy functional. Axially 



14l |. [15|, but they no longer correspond 



scheme, they also correspond to the global minima o 

symmetric hopfions A^i exist also for higher N = m 
.0 gl obal energy minima. For A^ = 3 the ground state hopfion is not planar and is called in 
158| A31, which can be viewed as deformed A31, with a pretzel-like position curve bent in 

3D to brake the axial symmetry. However, for A^ = 4 the axial symmetry is restored again 

in th e ground state, A22, which seems to have a similar to the A12 two-ring structure 90|, 



176| . The bent A41 also exists, but its energy it higher. 



Up to now all hopfions have been the simplest knots topologically equivalent to a circle, 
called unknots in the knot classification. A new phenomenon arises for A^ = 5, since the 
fundamental hopfion in this case consists of two linked unknots. This has nothing to do 
with the linking of preimages determining the value of A^. This time the position curve itself 
consists of two disjoint loops, corresponding to a charge 2 unknot linked to a charge 1 unknot. 
The Hopf charge is not simply the sum of charges of each component, but contains in addition 
the sum of their linking numbers due to their linking with the other components. It is worth 
noting that the linking number of the oriented circles can be positi ve or negative, depending 
on how they are linked 89|] (see FigJS]). Using the notation of Ref. 158l |. the A^ = 5 hopfion 





FIG. 5: Two possible ways to link A^ = 1 hopfions. Left: the total linking number is 2 = 1 + 1 and 
the total charge ifiV = l + l + 2 = 4. Right: the total linking number is —2 = — 1 — 1 and the 
total charge isA^ = 1 + 1 — 2 = 0. The numerical relaxation of these configurations gives therefore 
completely different results |59l ]. 



can be called ^{2, where the subscripts label the Hopf charges of the unknot components of 
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the link, and the superscript above each subscript counts the extra hnking number of that 
component. The total Hopf charge is the sum of the subscripts plus superscripts. Similarly, 
for A^ = 6 the ground state hopfion is a link of two charge 2 unknots, so that it can be called 

t1,1 
^2,2- 

For N = 7 the true knot appears at last: the ground state configuration corresponds in 
this case to the simplest non-trivial torus knot: trefoil knot. Let us remind that a {p, q) 
torus knot is formed by wrapping a circle around a torus p times in one direction and q times 
in the other, where p and q are coprime integers, p > q. One can explicitly parametrize 
it as p = -R + cos {p(f/q), z = sin {pip /q), where R > 1. A {p,q) torus knots can also be 
obtained as the intersection of the unit three sphere S"^ G C^ defined by 101^+ |crp = 1 with 
the complex algebraic curve o"^ + 0'^ = 0. The trefoil knot is the (3, 2) torus knot, and it 
determ ines the profile of the position curve of the N = 7 fundamental hopfion denoted 7K3 2 
in Ref. 



158 |. 



Sutcliffe (158| extends the energy minimization up to A^ = 16 specially looking for other 
knots. For the input configurations in his numerical procedure he uses fields parametrized 
by the rational map ansatz, 

W = V, 2.23 

where 0<aGZ, 0<6gZ and W,(f),a are defined by Eqs. fl2.15l) . fl2.16p . The position 
curve for such a field configuration coincides with the (p, q) torus knot, since the condition 

le k not equation. The parameters a, b determine the value of the Hopf 



W = 00 reproduces t 

charge, N = aq + bp 158l |. If p, q are not coprime, then the denominator in (12.231) factorizes 
and the whole expression describes a link. As a result, fixing a value of A^ one can construct 
many different knot or link configurations compatible with this value. Numerically relaxing 
these con figurations, Sutcliffe finds many new energy minima, discovering new knots and 



links 158| . He also obtains configurations that he calls x whose position curve seems to self- 
intersect and so it is not quite clear to what type they belong, unless the self-intersections 
are only apparent and can be resolved by increasing the resolution. 



The properties of all known hopfions, according to the results of Refs. 90|, 1581 ] . 15| . 
are summarized in Table [T] and in Figj6l Table [T] shows the Hopf charge, the type of the 
solution, with the ground state configuration in each topological sector underlined, and also 
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TABLE I: Known Hopf solitons according to Refs. 



^Q. 



ad], um 



N 


1 


2 


2 


3 


3 


4 


4 


4 


5 




Ai,i 


A2,l 


Ai,2 


A3,l 


A3,l 


A2,2 


A4,l 


A4,l 


t1,1 
^1,2 


£ 


1 


0.97 


0.98 


1.00 


1.01 


1.01 


1.03 


1.06 


1.02 


N 


5 


5 


6 


6 


6 


7 


7 


8 


8 




A5,l 


As,! 


t1,1 
^2,2 


t1,1 
^1,3 


Ae,! 


K3,2 


A7,l 


t1,1 
^3,3 


K3,2 


S 


1.06 


1.17 


1.01 


1.09 


1.22 


1.01 


1.20 


1.02 


1.02 


N 


8 


9 


9 


10 


10 


10 


11 


11 


11 




Ag.i 


T 2,2,2 
^1,1,1 


K3,2 


T 2,2,2 
^1,1,2 


t2,2 
^3,3 


K3,2 


T 2,2,2 
1,2,2 


K5,2 


t2,2 
^3,4 


£ 


1.40 


1.02 


1.02 


1.02 


1.02 


1.03 


1.02 


1.03 


1.04 


N 


11 


12 


12 


12 


12 


13 


13 


13 


13 




K3,2 


T 2,2,2 
■'-'2,2,2 


K4,3 


K5,2 


t2,2 
^4,4 


K4,3 


Xl3 


K5,2 


t3,3 
^3,4 


£ 


1.05 


1.01 


1.01 


1.04 


1.04 


1.00 


1.03 


1.04 


1.05 


N 


14 


14 


14 


15 


15 


15 


16 








K4,3 


K5,3 


K5,2 


Xl5 


T 4,4,4 
^1,1,1 


K5,3 


Xl6 






£ 


1.00 


1.01 


1.05 


1.01 


1.02 


1.02 


1.01 







the relative energy £ defined by the relation 



E^/E, = £N^'\ 



(2.24) 



where Ei is the energy of the A^ = 1 hopfion. Of the two decimal places of values of £ shown 
in the table the second one is rounded. Different groups give slightly different values for the 
energy, but one can expect the relative energy to be less sensitive to this. The values of £ for 
N< 7 shown in the table correspond to the data of Hietarinta and Salo 90| and of Sutcliffe 



158 



and it appears that for solutions described by both o 



the same. The data for 8 < A^ < 16 are given by Sutcliffe 



these groups these values are 



1581], apart from those for the 



Ajy,i hopfions for A^ = 5, 6, 7, 8, which are found in the earlier work of Battye and Sutcliffe 



151 ]. Although Atv,! solutions seem to exist also for A^ > 8, no data are currently available 



for this case. To obtain the energies from Eg. (12. 241) one can use the value Ei = 1.22, which 
is known to be accurate to the two decimal places 1721 ]. 

As one can see, the hopfion energies follow closely the topological lower bound (12.180 . 
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This suggests that the ground state hopfions actually attain this bound, so that they should 
be topologically stable. According to the data in Tabled one has inf{£^} = 0.97 for N < 16, 
and if this is true for all values of A^ then the optimal value for the constant in the bound 
(KWf is c = El inf{^} = 1.18. 



A rigorous existence proof for the hopfions was given by Lin and Yang 1191 ] , who demon- 
strate the existence of a smooth least energy configuration in every topological sector whose 
Hopf charge value belongs to an infinite (but unspecified) subset of Z. This shows tha t 
ground state hopfions exist, although perhaps not for any A^ G Z. As shown in Ref. 119| . 
their energy is bounded not only below but also above as 



E <C\N\ 



3/4 



(2.25) 



where C is an absolute constant. This implies that knotted solitons are energetically pre- 
ferred over widely separated unknotted multisoliton configurations when A^ is sufficiently 
large. Indeed, for a decay into charge one elementary hopfions the energy should grow at 
least as A^ for large A^, but it grows slower. 
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FIG. 6: Schematic profiles of the pos ition curves for the known hopfions (excepting the x-solutions) 



according to the results of Refs.[90],[l58|]. 



The position curves of the solutions, schematically shown in FigEl present an amazing 
variety of shapes. It should be stressed though that other characteristics of the solutions, as 
for example their energy density, do not necessarily show the same knotted pattern. Several 
types of links and knots appear, and each particular type can appear several times, for 
different values of the Hopf charge. Intuitively, one can view the position curves as wisps 
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made of two intertwined l ines corresponding to preimages of two infinitely close to — ric 



points on the target space 158| . Increasing the twist of the wisp increases the Hopf charge, 



without necessa rily changing the topology of the position curve. A more detailed inspection 



(see pictures in |l58l |) actually shows that configurations appearing several times in Fig|6l as 
for example the trefoil knot K3^2, become more and more distorted by the internal twist as A^ 
increases. Finally, for some critical value of A^, the excess of the intrinsic deformation makes 
it energetically favorable to change the knot/link type and pass to other, more complicated 
knot/link configurations. Estimating the length of the position curves shows that it grow s 



as A^^/^, so that the energy per unit length is approximately the same for all hopfions 158 |. 



E. Conformally invariant knots 

Hopf solitons in the Faddeev-Skyrme model, also sometimes called in the literature knot 
solitons of Faddeev-Niemi, of Faddeev-Skyrme, or of Faddeev-Hopf provide the best known 
example of knot solitons in field theory. However, there are also other field theory models 
admitting knotted solitons with a non-zero Hopf index. An interesting example proposed 



by Nicole 128l | is obtained by taking the first term in the Faddeev-Skyrme model (|2.ip and 



raising it to a fractional power, 

>CNicoie = {-d^n ■ d^nf^ . (2.26) 

A similar possibility, suggested by Aratyn, Ferreira and Zimerman (AFZ) |9(], uses the second 
term in the Faddeev-Skyrme Lagrangian, 

/:afz = iJ^^.J^'-'f' . (2.27) 

Both of these models are conformally invariant in three spatial dimensions and so the exis- 
tence of static solitons is not excluded for them by the Derrick argument. In fact, static knot 
solitons in these models exist and can even be obtained in close analytical form, which is 
achieved by simply using the axial ansatz (12.101) with the function 0, tp expressed in toroidal 
coordinates (12.131) according to Eq. (l2.14p |8|], |3[. Curiously, this separates away the v,(p 
variables in the field equations reducing the problem to an ordinary differential equation 
for Q{u) (in the Faddeev-Skyrme theory this does not work). In the AFZ model solutions 
for 6(m) can be expressed in terms of elementary functions for any n,m 8|], while in the 
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\7n\ 



1 case, where 



Nicole model they are obtained numerically [3j, apart from the \n\ 

;he s olution turns out to be the same in both models and is given by tan(0/2) = sinhti 

128| . 91]. These results remain, however, interesting mainly from the purely mathematical 

point of view at the time being, since it is difficult to justify physically the appearance of 

the fractional powers in Eqs.f l2.26ll . fl2.27l) . 

Other examples of solitons with a Hopf charge will be discussed below in Sec J VI C1 and 

Sec JVII C 21 



F. Can one gauge the knot solitons ? 



The Faddeev-Skyrme theory is a global field model, so that it cannot be a fundamental 
physical theory like gauge field theory models, but perhaps can be viewed as an effective 
theory. This suggests using the Faddeev-Skyrme knots for an effective description of some 
physical objects, and so it has been conjectured by Faddeev and Niemi that they could be 



used 



68| 



or a n effec tive description of glueballs in the strongly coupled Yang- Mills theory 66 ] 



153| . |l54| . [53|, [69|. This conjecture is very interesting, qui te in the spirit of the 



original Lord Kelvin's idea to view atoms as knotted ether tubes 162l |. and perhaps it could 
apply in some form. In fact, when describing the ?7(1440) meson, the Particle Data Group 
says (see p. 591 in Ref . J60| ) that "the ?7(1440) is an excellent candidate for the 0""*" glueball 
in the flux tube model [691]". 

Some other physical applications of the global fleld theory knot solitons could perhaps 
be found. However, if they could be promoted to gauge fleld theory solutions, then they 
would be much more interesting physically, since in this case they would flnd many inter- 
esting applications, as for example in the theories of superconductivity and of Bose-Einstein 
cond ensation [ill], [l2|, in the theory of plasma [67|, [631], in Standard Model [361], |70| . 
130| . o r per haps even in cosmology, where they could presumably describe knotted cosmic 
strings 1651 ]. For this reason it has been repeatedly conjectured in the literature that some 
analogs of the Faddeev-Skyrme knot solitons could also exist as static solutions of the gauge 
fleld theory equations of motion. This conjecture is essentially inspired by the fact that the 
Faddeev-Skyrme theory already contains something like a gauge fleld: JF^^,. Moreover, we 
shall now see that changing the variables one can rewrite the theory in such a form that it 
looks almost identical to a gauge fleld theory (or the other way round). 
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G. Faddeev-Skyrme model versus semilocal Abelian Higgs model 

Let $ be a doublet of complex scalar fields satisfying a constraint, 

$= I '^ I , $t$ = |0|2 + |^|2 = i^ (2.28) 

such that ^ & S^. Let us consider a field theory defined by the Lagrangian density 

£[$] = -J J'^.T^"' + (2)^$)tl?^$ (2.29) 

with jF^i, = dfj,Au — dyA^ and P^$ = {d^ — iA^)^ and with 

A^ = -i<^^d^^. (2.30) 

In fact, this theory is again the Faddeev-Skyrme model but rewritten in different variables, 
since upon the identification 

^a ^ ^]^a^ ^2.31) 



(r° being the Pauli matrices) the fields A^j, and jF^j/ coincide with those in (12. ip 14 1| and 



the whole action (I2.29P reduces (up to an overall factor) to (12. ip 65l |. More precisely, (I2.3ip 
is the Hopf projection from S^ parametrized by (0, a) to S'^ parametrized by the complex 
projective coordinate 0/cr. For example, the axially symmetric fields (I2.10p . (12. lip are 
obtained in this way by choosing the CP^ variables 

(^ = cos - e^'^'^, (J = sin - e^"*^, (2.32) 

such that the phases of </> and a wind, respectively, along the two orthogonal direction as 
shown in Figl2] and the Hopf charge is A^ = nm. 

The fields in the static limit, $ = $(x), can now be viewed as maps S'^ -^ S^ , but their 
energy 

E[<^] = j (\Vk^^ + ^ {T.,kf\ t/^x (2.33) 

is still bounded from below as in Eq. (l2.18p . The topological charge A^ = A^[$], still expressed 
by Eq. (l2.9p . is now interpreted as the index of map S^ ^ S'^ , N E 7r3(5'^) = Z. The theory 
therefore admits the same knot solitons as the original Faddeev-Skyrme model. However, in 
the new parametrization the theory looks almost like a gauge field theory, in particular it 
exhibits a local U(l) gauge invariance under $ -^ e*"$, A^ ^> A^ + d^a. 
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Let us now compare the model fl2.29p to a genuine gauge field theory with the Lagrangian 
density 

£[$, A,] = -i F^^F'^'' + {D^^^D^^ - ^ (<l>t$ - l)^ . (2.34) 

Here $ is again a doublet of complex scalar fields, but this time without the normalization 
condition (12.281) . the condition (12.301) being also relaxed, so that A^ is now an independent 
field. One has F^^ = d^A^ — d^A^ and D^$ = {d^ — iA^)^. This semilocal [2| Abehan 
Higgs model with the SU(2)giobaixU(l)iocai internal symmetry arises in different contexts, 
in particular it can be viewed as the Weinberg-Salam model in the limit where the weak 
mixing angle is 7r/2 and the SU(2) gauge field decouples. The non-relativistic limit of this 



theory is the two-component Ginzburg-Landau model |77l |. 

Let us now consider the limit A ^ oo. In this sigma model limit the constraint (I2.28P 
is enforced and the potential term in (I2.34p vanishes. The theories (I2.29P and (I2.34p then 
look identically the same, the only difference being that in the first case the vector field Afj, 
is defined by Eq. (l2.30p and so is composite, while in the second case A^ in an independent 
field. 

H. Energy bound in the Abelian Higgs model 

The question now arises: does the gauge field theory (I2.34p . at least in the limit where 
$^$ = 1, admit knot solitons analogues to those of the global model (I2.29P ? If exist, such 
solutions would correspond to minima of the energy in the theory (12.340 in static, purely 
magnetic sector, 

E[$, Ak] = J (\Dk<^\' + ^ {Fik)A rf'x . (2.35) 

At first thought, one may think that the answer to this question should be affirmative. 
Indeed, the energy functionals (I2.33P and (12.350 look identical. They have the same internal 
symmetries and the same scaling behaviour under x -^ Ax. The two theories also have the 
same topology associated to the field $, since in both cases $(x) defines a mapping S^ -^ S^ 
with the topological index (12.90 . 

The gauged model (12.350 contains in fact even more charges than the global theory (12.330 . 
since it actually has two vector fields: the independent gauge field A^ and the composite 
field Ak = ^{dk^^^ — $'''(9^$). It is convenient to introduce their difference C^ = A^ — Ak- 
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Defining the linking number between two vector fields, 

I[A, B] = 2^f e,,kAd,Bk d'^ , (2.36) 

one can construct three different charges, 

iV[$] = I[A,A], L = I[C,A], Ncs[A] = I[A,A], (2.37) 

which are, respectively, the topological charge (12. 9p . the linking number between A^ and Ck, 
and the Chern-Simons number of the gauge field. The following inequality, established by 



Protogenov and Verbus 



E[$,A,]>c,|iV|3/^(^l-Hy, (2.38) 



holds: 

2 



where ci is a positive constant. This can be considered as the generalization of the topological 
bound (12.181) of Vakulenko-Kapitanski. 

Given all these, one may believe that the local theory (12.351) admits topologically stable 
knot solitons similar to those of the global theory (12.331) . 

I. The issue of charge fixing 

The difficulty with implementing the Protogenov- Verbus bound (I2.38P in practice is that 
it contains two different charges, N and L, but without invoking additional physical assump- 
tions there is no reason why L should be fixed while minimizing the energy. 

Let us consider first the charge A^ = A^[$]. Its variation vanishes identically, so that 
it does not change under smooth deformations of $. It is therefore a genuine topological 
charge whose value is completely determined by the boundary conditions, and so it should 
be fixed when minimizing the energy. 

Let us now consider the linking number L=/[C, A]. The analogs of this quanti ty ha ve been 
studied in the theory of fluids, where they are known to be integrals of motion 178| . In the 
context of gauge field theory, this quantity is gauge invariant. However, it is not a topological 
invariant, since its variation does not vanish identically and so it does change under arbitrary 
smooth deformations of Ck = Ak — Ak- On cannot fix L using only continuity arguments, 
because there are no topological conditions imposed on Ak-, whose arbitrary deformations 
are allowed. The only topological quantity associated to Ak could be a magnetic charge 
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related to a non-trivial U(l) bundle structure. However, since we are interested in globally 
regular solutions, the bundle base space is M^ (or S^) without removed points, in which case 
the bundle is trivial. 

As a result, on continuity grounds only, one can fix A^ but not L. But then, as is obvious 
from Eg. (12.381) . there is no non-trivial lower bound for the energy, since one can always 
choose L=A^ in which case the expression on the right in (I2.38P vanishes. More precisely, 
since there are no constraints for A^, nothing prevents one from smoothly deforming it to 
zero, after which one can scale away the rest of the configuration. Explicitly, given fields 
$, Ak one can reduce -E[$, Ak] to zero via a continues sequence of smooth field deformations 
preserving the value of the topological charge N[^], 

E[<D(x), Afc(x)] -^ E[<l>(Ax),7Afc(x)] (2.39) 

by taking first the limit 7 — > and then A — > cxo J74l ]. 

The conclusion is that without constraining Ak, with only A^[$] fixed, the absolute min- 
imum of £'[$,y4fc] is zero, so that there can be no absolutely stable knots. To have such 
solutions, one would need to constraint somehow the vector field A^, for example it would 
be enough to insure that Cu = Au — Ak be zero or small. Such a condition is often assumed 

n n h [^ 

in the literature [12], [ll|, [36[, [70|], but usually simply ad hoc. Unfortunately, it cannot be 
justified on continuity grounds only, without an additional physical input. 



J. Searching for gauged knots 



The above arguments do not rule out all solutions in the theory. Even though the global 
minimum of the energy is zero, there could still be non-trivial local minima or saddle points. 
The corresponding static solutions would be metastable or unstable. One can therefore 
wonder whether such solutions exist. This question has actually a long history, being first 



addressed by de Vega [49| and by Huang and 
repeatedly reconsidered by different authors 



ipton |93[| over 30 yea rs ago , and being then 



in4 |i3nt . y, [IT^, M, [stI, y, |56| 



However, the answer is still unknown. No solutions have been found up to now, neither has 
it been shown that they do not exist. 

Trying to find the answer, all the authors were minimizing the energy given by the sum 
of -E'[$,ylfc] and of a potential term that can be either of the form contained in (I2.34p or 
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a more general one. A theory with two ve ctor fields with a local U(l)xU(l) invariance 



Witten's model of superconducting strings 177 1 - has also been considered in this context 

n n n 

|57| . [551], [56!]. The energy was minimized within classes of fields with a given topological 
charge A^ and having profiles of a loop of radius R. The resulting minimal value of energy 
was always found to be a monotonously growing function of -R, thus always showing the 
tendency of the loop to shrink thereby reducing its energy. 

These results render the existence of solutions somewhat implausible. However, they do 
not yet prove their absence. Indeed, local energy minima may be difficult to detect via energy 
minimization, as this would require starting the numerical iterations in their close vicinity, 
because otherwise the minimization procedure converges to the trivial global minimum. In 
other words, one has to choose a good initial configuration. However, since 'there is a lot of 
room in function space', chances to make the right choice are not high. 



It should be mentioned that a positive result was once reported in Ref. |l30l ]. were the 
energy was minimized in the N = 2 sector and an indication of a convergence to a non-trivial 
minim um w as observed. However, this result was not confirmed in an independent analysis 



in Ref. 174| . so that it is unclear whether it should be attributed to a lucky choice of the 



initial configuration or to some numerical artifacts. 

III. KNOT SOLITONS IN GAUGE FIELD THEORY 

As we have seen, fixing only the topological charge does not guarantee the existence of 
knot solitons in gauge field theory. In order to obtain such solutions one needs to constraint 
the gauge field in order that it could not be deformed to zero. Below we describe two known 
examples of such solutions. 

A. Anomalous solitons 

One possibility to constraint the gauge field is to fix its Chern-Simons number. An 



exam ple of how this can be done was suggested long ago by Rubakov and Tavkhelidze 143| , 



142| . who showed that the Chern-Simons number can be fixed by including fermions into 
the system. They considered the Abelian Higgs model (12.341) . but with a singlet and not 
doublet Higgs field, augmented by including chiral fermions. In the weak coupling limit 
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at zero temperature this theory contains states with A'"^ non-interacting fermions and with 
the bosonic fields being in vacuum, A^ = 0, $ = 1. The energy of such states is -E ~ Np. 
Rubakov and Tavkhehdze argued that this energy could be decreased via exciting the bosonic 
fields in the following way. 

Owing to the axial anomaly, when the gauge field A^ varies, the fermion energy levels 
can cross zero and dive into the Dirac see. The fermion number can therefore change, but 
the difference Np — Nq^ is conserved. As a result, starting from a purely fermionic state 
and increasing the gauge field, one can smoothly deform this state to a purely bosonic state, 
whose Chern-Simon number will be fixed by the initial conditions. 



iVcs = Np. (3.40) 



Now, the energy of this state. 



E[^, A,] = J (^|D,$r + \ {F,uf + \m'- 1)') d'^ , (3.41) 

can be shown to be bounded from above by i?o(^cs)'^^^ where Eq is a constant, and so for 
arge Np = Ncs it grows slower than the energy of the original fermionic state, E ~ A'^cs 
143| . 142| . Therefore, for large enough Np, it is energetically favorable for the original 



yurely fermionic state to turn into a purely bosonic state. The latter is called anomalous 



143| , 142l | . The energy of this anomalous state can be obtained by minimizing the functional 
(I3.4ip with the Chern-Simons number fixed by the condition (I3.40p . 

Suc h an energy minimization was carried out in the recent work of Schmid and Shaposh- 



nikov 



1511 ]. First of all, they established the following inequality, 

E[<^,Ak]>c{Ncs?/\ (3.42) 



which reminds very much of the Vakulenko-Kapitanski bound ( J2.18P for the Faddeev-Skyrme 
model, but with the topological charge replaced by A'cs- This gives a very good example of 
how constraining the gauge field can stabilize the system: even though in the theory with a 
singlet Higgs field there is no topological charge similar to the Hopf charge, its role can be 
taken over by the Chern-Simons charge. 

In order to numerically minimize the energy, Schmid and Shaposhnikov considered the 
Euler- Lagrange equations for the functional 

E[^, Ak] +1^1 e.jkAd.Akd^x (3.43) 
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FIG. 7: Spindle torus shape of the anomalous solitons for large A'^cs- 

where £"[$, Ak] is given by (13.411) and fi is the Lagrange multiplier. In the gauge where 
$ = $* = these equations read (with A = A^) 



A(j)-A- 



\i^^ 



V X (V X A) + 2/iV X A + 2(j)^A = 0, 



(3.44a) 
(3.44b) 



where V and A = (V)^ are the standard gradient and Laplace operators, respectively. 
Multiplying Eq. (I3.44bp by A and integrating gives the expression for the Lagrange multiplier, 

^ ' ' '- ■ ■■" (3.45) 



/^ 



327r2iVcs 



(V X AY + \A\'(j)' d-'x 



Solutions of the elliptic equations (I3.44ap . (13.44bl) were studied in |l5l| in the axially sym- 
metric sector, where = 0(p, z) and A = A{p, z), with the boundary conditions 



at infinity and 



1, A = 



dp<j) = dpA, = 0, Ap = A^ 







(3.46) 



(3.47) 



at the symmetry axis p = 0. The solutions obtained are quite interesting. They are very 
strongly localized in a compact region, Q, of the (p, z) plane centered around a point (p^, 0). 
Inside Q the field A^ is non-zero, while is almost constant and very close to zero. As 
one approaches the boundary of the region, dQ, the field Aj^ tends to zero, while is still 
almost zero. Finally, in a small neighborhood of dQ whose thickness is of the order of the 
Higgs boson wavelength, starts varying and quickly increases up to its asymptotic value. 
Outside fi one has everywhere 0^1 and A^ ^ 0, such that the energy density is almost 
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zero. The energy for these solutions scales as (A^cs)^ • For large A'^cs the region fl can be 
described by the simple analytic formula, 

{p-psf + z^KRl (3.48) 

where Rs > Ps- In other words, the 3D domain where the soliton energy is concentrated can 
be obtained by rotating around the 2;-axis a disc centered at a point whose distance fro m th e 



z-axis is less than its radius (see Fig|7]). Such a geometric figure is called spindle torus 151 |. 
The solutions for large A^'cs are very well approximated by setting = 1 and A^ = outside 
the spindle torus, while inside it one has = and Ak is obtained by solving the linear 
equation (I3.44bl) . The spindle torus approximation becomes better and better for large A'^cs, 
in which limit Schmid and Shaposhnikov obtain the following asymptotic formulas for the 
energy and parameters of the torus, which agree very well with their numerics, 

1 7 /TV \ ^^^ 
E = 118X'/'Mw{Ncsf\ p^ = __(^j , i?, = 1.49p,, (3.49) 

where Mw is the vector field mass. 

It is likely that these solution attain the lower energy bound (I3.40p . which means that 
they should be topologically stable. It should, however, be emphasized that the anomalous 
solitons require a rather exotic physical environment, since for them to be energetically 
favoured as compared to the free fermion condensate, the density of the latter should attain 
enormous values possible perhaps only in the core of neutron stars. 

Summarizing, fixing the Chern-Simons number forbids deforming the gauge field to zero 
and gives rise to stable knot solitons in gauge field theory, even if the Higgs field is topolog- 
ically trivial. It is unclear whether this result can be generalized within the context of the 
model (I2.34P with two component Higgs field - since the Protogenov-Verbus formula (I2.38P 
contains not the Chern-Simons number but the linking number L. 

B. Non-Abelian rings 

Another interesting class of objects arises in the non-Abelian gauge field theory, where 
one can have smooth, finite energy loops stabilized by the magnetic energy. 
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1. Yang- Mills- Higgs theory 

Let us parametrize the non-Abelian Yang-Mills-Higgs theory for a compact and simple 
gauge group Q as 

£[A^,$] = --^{F^,Fn + {D^^)^D^'^ - f/($) . (3.50) 

Here the gauge field strength is F^^ = d^Ay - d^A^ - ig[Af,, A^] = F^JTa where A^ = A^T^ 
is the gauge field, a = 1, 2, ... , dim(^), and g is the gauge coupling. The Hermitian gauge 
group generators T^ satisfy the relations 

[T„, Tb] = iCabcTc, tr(T,Tb) = KSab ■ (3.51) 

The Lie algebra inner product is defined as {AB) = ^tT{AB) = A°'B"'. The Higgs field $ 
is a vector in the representation space of Q where the generators T^ act; this space can be 
complex or real. The covariant derivative of the Higgs field is -0^$ = (9^j — igA^^)^ and the 
Higgs field potential can be chosen as t/($) = ^ ($'''$ — 1)^. The Lagrangian is invariant 
under the local gauge transformations, 

$ ^ U<l>, A^ -. V{A^ + -d^)V-\ (3.52) 

where U = exp{ia°-{x^^)Ta) G Q. The field equations read 



V^Ff"" = ig {{D^^yTa^ - ^^TaD"^} T„ , 

dU 
9(¥ti>) 



D^DM$ = ___--_<|., (3.53) 



where P^ = d — ig[Afj_, ] is the covariant derivative in the adjoint representation. The 
energy momentum tensor is 

T^' = -I^Fy^F^"") + (D^$)tD'^$ + {D^'^^D^^ - 5>^C . (3.54) 

2. Monopole rings 



Th e rin g soli tons in the theory fl3.50p have been first constructed by Kleihaus, Kunz and 



Shnir 



106(1, 107l | in the case where ^=SU(2) and the Higgs field is in its adjoint, $ = $", such 



that the gauge group generators are 3x3 matrices with components (Ta)bc = —i^ahc- The 
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und a men tal solutions in this theory are the magnetic monopoles of 't Hooft and Polyakov 



160| . 136l |. while the ring solitons are more general solutions. Specifically, in the static, 



axially symmetric and purely magnetic case it is consistent to choose the following ansatz 
for the fields in spherical coordinates: 

A^dx^ = {Kidr + (1 - K2)d^)T^ + m{K3Tr + (1 - K^)Ts) sin^ dip, 
T„$" = 0iT, + 02T^ , (3.55) 



where functions Ki, K2, K3, K4, 0i, 02 depend o n r, -i ^ anc 
conditions at the symmetry axis and at infinity 106l |. 107l | 



are subject of suitable boundary 
. Here 



sin{k'd) cos{rrnp)Ti + sin(A;'i9) sm{mip)T2 + cos{ki})T3 
1 d 



T,9 — — -— -Ti- 






k d'd ^ '' "^ m sin •§ dip 

with k,m E Z. Using the gauge invariant tensor 



(3.56) 



flU 



/lU 



(3.57) 



and its dual, JF^,^ = |e^,^a/3-^a/3) one can define the electric and magnetic currents, respec- 
tively, as 

J>. = d-^Ta^, j^ = d''Taf.. (3.58) 

It turns out that the solutions depend crucially on values of k,m in (13.551) . In particular, 
their magnetic charge is given by 



^ ITT- r 



(3.59) 



The following solutions are known in the limi t of v anis hing Higgs potential (for a generic 

potential their structure is more complicated) 1061 ]. 107|: 

k = l,m = 1 - the spherically symmetric 't Hooft-Polyakov monopole. 

k = l,m > 1 - multimonopoles. 

k > l,m = 1,2 - monopole-antimonopole sequences. 

k = 2l>2,'m>3- monopole vortex rings. 

In the first three cases the Higgs field has discrete zeros located at the z-axis. Solutions of the 

last type are especially interesting in the context of our discussion, since zeros of the Higgs 

fields in this case are not discrete but continuously distributed along a circle (for k = 2) 
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FIG. 8: Left: Higgs field amplitude for the k = 2, m = 3 monopole ring solution in the limit 
where the Higgs field potential is zero. \^\ vanishes at a point in the {p,z) plane away from the 
z-axis, which corresponds to a ring. Right: schematic shape of charge/current distribution for this 
solution. 

around the 2;-axis (see FigJHD- Solutions in this case can be visualized as stationary rings 



stabilized by the magnetic energ y. T 
and can be elucidated as follows 155 



l e me chanism of their stabilization is quite interesting 



1561 ]. If one studies the profiles of the currents ()3.58p 



for these solutions, it turns out that both the magnetic charge density jo and the electric 
current density jk have ring shape distributions, as qualitatively shown in FigJHl 

Although the total magnetic charge is zero, locally the charge density is non-vanishing 
and the system can be visualized as a pair of magnetically charged rings with opposite charge 
located at z = ±zq, accompanied by a circular electric current in the z = plane. The two 
magnetic rings create a magnetic field orthogonal to the z = plane. This magnetic field 
forces the electric charges in the plane to Larmore orbit, which creates a circular current. 
The Biot-Savart magnetic field produced by this current acts, in its turn, on the magnetic 
rings keepin g th em away from each other, so that the whole system is in a self-consistent 
equilibrium 155| (assuming the magnetic rings to be rigid). 

It is, however, unlikely that this sophisticated balance mechanism stabilizing the rings 



against contraction could also guarantee their stability with respect to all possib 



tions. In fact, the monopole-antimonopole solution is known to be unstable 161| . while the 



e deforma- 



monopole rings can be viewed as generalizations of this solution. They are therefore likely 
to be saddle points of the energy functional and so they should be unstable as well. 
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3. Sphaleron rings 
Very recently, a similar ring construction was carried out by Kleihaus, Kunz and Leissner 



101| within the context of the Yang-Mills-Higgs theory (13.501) with ^=SU(2) and with the 
Higgs field in its fundamental complex doublet representation, where Tq = ^Tq. This theory 
can be viewed as the SU(2) xU(l) Weinberg-Salam theory in the limit where the weak mixing 
angle vanishes and the U(l) gauge field decouples. Kleihaus, Kunz and Leissner used exactly 
the same field ansatz fl3.55p . only modifying the Higgs field as 

$ = (0iT, + 02T^)$o , (3.60) 

where $o is a constant 2-vector. As in the monopole case, in this case too the solutions 



depend strongly on the choice of the integers k and m in Eg. (13.561) . The fundamental 
solutions in this case are the sphalerons - unstable saddle point configurations that can be 
smoothly deformed to vacuum. They can be characterized by the Chern-Simons number, 
given by the same formula as the magnetic charge in the monopole case, up to the factor 
1/2, so that half-integer values are now allowed: Q = m[l — ( —1)^ 1/4. 



Setting k = m = 1 gives the Klinkhamer-Manton sphaleron |l08l |. in which case the Higgs 
field vanishes at one point. Choosing /c = l,?7i>lor/i;>l,?72 = l gives multisphalerons or 
sphaleron-antisphalero n so lutions for which the Higgs field has several isolated zeros located 



at the symmetry axis [lOl|. A new type of solution arises for /c > 2, ?7i > 3, in which case 
the Higgs field vanishes on one or more rings centered around the symmetry axis. In this 
respect these sphaleron rings are quite analogues to the monopole rings. It is unclear at the 
moment whether their existence can be qualitatively explained by a mechanism similar to 
that for the monopole rings, shown in FigJHl 

Since sphalerons are unstable objects, it is very likely that sphaleron rings are also un- 
stable. However, similar to the sphalerons, they could perhaps be interesting physically as 



mediators of baryon number violating processes 108 1. 



IV. ANGULAR MOMENTUM AND RADIATION IN FIELD SYSTEMS 

We are now passing to the spinning systems with the ultimate intention to discuss spin- 
ning vortex loops stabilized by the centrifugal force - vortons. As was already said in the 
Introduction, usually vortons are considered within a qualitative macroscopic description as 
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loops made of vortices and stabilized by rotation. This description is suggestive, but it does 
not take into account the radiation damping. At the same time, the presence of the vorton 
angular momentum requires some internal motions in the system (see Figj9]), as for example 
circular currents, and these are likely to generate radiation carrying away both the energy 
and angular momentum. It is therefore plausible that macroscopically constructed vortex 
loops will not be stationary field theory objects, but at best only quasistationary, with a 
finite lifetime determined by the radiation rate. 

RADIATION ? 





vorton 



vorton 



FIG. 9: One can make a loop from a vortex carrying a current j and momentum P. It will have 
an angular momentum, but it may be radiating. 

The best way to decide whether vortons are truly stationary or only quasistationary is 
to explicitly resolve the corresponding field theory equations. The current situation in this 
direction is not, however, very suggesti ve. W ithin the original local U(l)xU(l) Witten's 
model of superconducting cosmic strings 177l | vorton solutions have never been constructed. 



edge , the only explicit vorton solutions have been presented by 



1181 ] within the global version of Witten's model, and also by 



To the best of our know 

Lemperier and Shellard 

Battye, Cooper and Sutcliffe [17[ in a special limit of the same global model. In addition, 

vortons in a 2 + 1 dimensio nal fi eld theory toy model have been recently analyzed 19|. 

Lemperier and Shellard J118l | considered the full hyperbolic evolution problem for the 
fields, with the initial data corresponding to a vortex loop. Evolving dynamically this loop 
in time, they saw it oscillate around a visibly stationary equilibrium position, and they could 
follow these oscillations for several dozens characteristic periods of the system. This suggests 
that vortons exist and are stable against perturbations. However, one cannot decide on these 
grounds whether the equilibrium configurations are truly stationary or only quasistationary, 
since the radiation damping could be non-zero but too small to be visible in their numerics. 
In addition, Lemperier and Shellard did not actually consider precisely the global version 
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of Witten's model (see Eg. (16.1 38 p below), but, in order to improve the numerics, added to 
it a Q-ball type interaction term of the form |(;/)|^|crp, where 0,0" are the two scalars in the 
model. 

Battye, Cooper and Sutcliffe [17| did not study Witten's model but minimized the energy 
of a non-relativistic Bose-Einstein condensate, which seems to be mathematically equivalent 
to solving equations of Witten's model in a special limit. They found non-trivial energy 
minima saturated by configurations of vorton type, which again suggests that vortons exist, 
at least in this limit. Moreover, this suggests that they are indeed non-radiative - since 
being already in the energy minimum they cannot loose energy anymore. It would therefore 
be interesting to construct these solutions in a different way, extending the analysis to the 
full Witten's model. 

The method we shall employ below to study truly non-radiating vortons will be to con- 
struct them as stationary solutions of the elliptic boundary value problem obtained by 
separating the time variable. However, first of all we need to understand how in principle 
a non-radiating field system can have a non-zero angular momentum. Both angular mo- 
mentum and radiation are associated to some internal motions in the system, and it is not 
completely clear how to reconcile the presence of the former with the absence of the latter. 

A. Angular momentum for stationary solitons 

In what follows we shall be considering field theory systems obeying the following four 
conditions: 

(1) stationarity 

(2) finiteness of energy 

(3) axial symmetry 

(4) non-vanishing angular momentum 

The angular momentum is defined as the Noether charge associated to the global space- 
time symmetry generated by the axial Killing vector K = d/d(f, 

J= fr^d^x. (4.61) 

Let us discuss the first three conditions. 

(1) A system is stationary if its energy momentum tensor Tj^ does not depend on time. 
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According to the standard definition of symmetric fields [7l|, for stationary fields the action 
of time translations can be compensated by internal symmetry transformations. 

If all internal symmetries of the theory are local, then there is a gauge where the compen- 
sating symmetry transformation is trivial, so that the stationary fields are time- independent. 
We shall call them manifestly stationary. If the theory contain also global internal symme- 
tries and if the compensating symmetry transformation is global, then its action cannot be 
trivialized and so the action of time translations will be non-trivial. The fields will explicitly 
depend on time in this case, typically via time- dependent phases, and we shall call them 
non-manifestly stationary. 

For example, in a system with two complex scalars coupled to a U(l) gauge field one 
cannot gauge away simultaneously phases of both scalars. A non-manifestly stationary field 
configuration will then be 0i(x), </)2(x)e*'^*, yl^(x). However, if there is only one scalar, then 
it is always possible to gauge away its time-dependent phase. 

(2) Even if Tjf is time- independent, one can still have a constant radiation flow compen- 
sated by the energy inflow from infinity. However, if the energy is finite, then the fields 
fall-off fast enough at infinity to eliminate this possibility. 

In principle, one can also have situations where T^^ is tir ne-d ependent, but radiation is 



13911. However, such cases are 



nevertheless absent, as for the breathers in 1-1-1 dimensions 
probably less typical and we shall not discuss them. 

(3) It is intuitively clear that asymmetric spinning systems will more likely radiate than 
symmetric ones. It is therefore most natural to assume spinning non-radiating solitons to 
be axially symmetric. The axially symmetry can be manifest or non-manifest. In fact, it 
appears that the axial symmetry condition can sometimes be relaxed, but such a possibility 
seems to be more exotic and will be discussed below only very briefly in Sec JVCi 

Let us now analyze possibilities for condition (4) to coexist with (l)-(3). 

1. The case of manifest symmetries - no go results 

In theories where all internal symmetries are local the stationary field are time- 
independent. Can one have an angular momentum in this case ? 

In fact, even without an explicit time dependence one can have a non- vanishing field 
momentum expressed by the Poynting vector, S x B, and this could give a contribution to 
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the angular momentum, 

fx{SxB)d^x. (4.62) 

If the theory contains gauged scalars, they will give an additional contribution. If there exists 
a stationary, globally regular, finite energy on-shell configuration for which the integral fl4.6ip 
is non-zero, this would correspond to a rotating and non-radiating soliton. The existence of 
such solutions is not a priori excluded. However, it turns out that for a number of physically 
interesting cases they can be ruled out. 

Specifically, since the time translations and spatial rotations commute, the system is not 



only manifestly stationary but also manifestly axially symmetric 71[. It turns out that 
the latter condition allows one to transform the volume integral in Eqs. fl4.6ip to a surface 
integral, which is often enough to conclude that J = 0. 

To me more precise, let us consider the Yang-Mills-Higgs theory (13.501) . If the fields 
(y4^, $) are manifestly invariant under the action of a Killing symmetry generator K = d/ds, 
then there is a gauge where they do not depend on the corresponding spacetime coordinate 
s and their Lie derivatives along K vanish, CkA^^ = Cj^^ = 0. In some other gauge the 



fields could depend on s, but only in such a way that 71 1 



CKA^ = -D^WiK), CK<^ = igWiK)^. (4.63) 

Here W{K) takes its values in the Lie algebra of the gauge group and transforms as connec- 
tion under gauge transformations, li K = d/d(f is the axial Killing vector and W^p = W{K), 
then in spherical coordinates these conditions reduce to 

d^A^ = V^W^ , 9^$ = igW^^ . (4.64) 

If these conditions are fulfilled, then it is straightforward to check with the field equations 
(I3.53P that the T? component of the energy- momentum tensor (I3.54p has a total derivative 
structure [50|, [l69|, 

Tl = du{{Wp - A^)F'^) . (4.65) 

Since both A^p and W^p transform as connections under gauge transformations, their differ- 
ence is gauge covariant, so that this formula is gauge invariant. For globally regular solutions 
the electric field £^ = F^'' is everywhere bounded, which allows one to transform the volume 
integral of T° to a surface integral over the boundary of M"^ at infinity. This gives the surface 
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integral representation for the angular momentum 



J 



{{W^-A^)S)dS. 



(4.66) 



This formula imposes rather strong restrictions on the existence of spinning solitons. For 
example, it shows that if the electric field S decays at infinity faster than 1/r^, such that 
the electric charge is zero, then J = 0. A non-zero electric charge is therefore necessary to 
have a non-zero angular momentum 50||. Similarly, using the asymptotic conditions in the 
far field zone, one can show that for a number of important cases the asymptotic behavior 
of the fields does not allow for the surface integral to be non-zero. 




FIG. 10: Do static solitons admit stationary, spinning generalizations ? 



Specifically, one can consider any known static, spherically symmetric soliton, as for 
example the 't Hooft-Polyakov monopole, and ask if it admits stationary, spinning general- 
izations (see Fig JTOl) ? For example, the static Schwarzschild black hole solution of Einstein's 
equation in General Relativity can be generalized to the manifestly stationary, spinning Kerr 
black hole. One can wonder if similar generalizations are possible for solitons in non-linear 
field theories in Minkowski space. If they exist, then in the far field zone their fields should 
approach those of the original static soliton, so that they should be spherically symmetric, 
up to small corrections which could contribute to the surface integral (I4.66p . Since the 
correction s are small, their most general form can be determined from the linearized field 
equations |l69 |. 

Surprisingly, one discovers in this way that none of the well known solitons in the gauge 
field theory ( 13.50p wit h gau ge group SU(2), as for example the magnetic monopoles of 
't Hooft-Polyakov 160||, 1361 ]. dyons of Julia-Zee 98|, and sphalerons of Klinkhamer-Manton 
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admit spinninff generalizations within the manifestly stationary and manifestly ax- 



isymmetric sector 50| , 5l| , 169| . The picture in Fig{TO] thus does not apply for these 
solitons. If the gauge symmetry is completely broken in the Higgs vacuum, as for the 
sphalerons, then the fields approach their vacuum values exponentially fast and the surface 
integral vanishes. For the monopoles and dyons there is a long-range field associated with 
the unbroken U(l), so that some additional analysis is required to show that in the asymp- 
totic reg ion t here are actually no field modes giving a non-zero contribution to the surface 



integral [169 1. 



As a result, the best known SU(2) solitons cannot spin in the manifestly stationary and 
manifestly axisymmetric sector. This restricts rather strongly the existence of spinning soli- 
tons with manifest symmetries, although does not rule them out completely. Such solutions 
might exist in theories with other gauge groups. Their explicit examples in U(l) gauge field 
theories will be presented below. Non-manifest ly symmetric rotational excitations could 
perhaps exist for the sphalerons, since only an SU(2) part of their U(2) internal symmetry is 
gauged. A more exotic possibility to have stationary rotation without axial symmetry will 
be briefly discussed in Sec JVCI 

2. The case of non-manifest symmetries 

More general possibilities to have spinning solitons arise in theories where not all internal 
symmetries are local, since in this case the conditions fl4.63p can be generalized as 

CkA^ = V^W{K) , Ck^ = tgW{K)^ + iT{K)^ , (4.67) 

where T(K) is a function of the global symmetry generators. One can gauge away W{K) 
but not T{K), so that the action of the Killing symmetry will always be non-trivial, since 
$ will depend explicitly on the corresponding spacetime coordinate s as 

$(s) =e'^^(^)$o. (4.68) 

This means that the s-dependence is equivalent to a sequence of internal symmetry trans- 
formations, in which case the invariant objects like T^ will not depend on s at all. Such 
a symmetry can be called non-manifest. Non-manifest ly stationary and non-manifestly ax- 
isymmetric fields would typically depend on t, ip via the complex phase factor 

exp{i{ujt -\- rmp)} , (4.69) 
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in which case one can say that it is the phase that spins. The angular momentum in this 
case cannot be totally expressed by a surface integral and will contain a volume integral 
contribution, in which case typically J ~ ujm. 

As a result, one can have J 7^ 0. However, the absence of radiation is not yet guaranteed. 
Separating the time variable makes the equations elliptic, but with the mass term(s) modified 
as 

M^ -^ M^ - tu^, (4.70) 

where M^ collectively denotes masses of the field excitations in the asymptotic zone. Now, 
if the equations admit globally regular solutions with 

0^2 < M^ (4.71) 



then these solutions will behave asymptotically as exp{ — vM^--u^ r} and there will be 
no radiation. On the other hand, if uo'^ > M^ then solutions in the asymptotic region will 
oscillate as exp{±.i\/uj'^ — M"^ r} thus showing the presence of the ingoing and outgoing 
radiation, even though T^ is time-independent, so that the total energy will be infinite. As 
a result, the no-radiation condition (14.711) becomes crucial in this case. 

3. Spinning solitons as solutions of elliptic equations 

Summarizing the above discussion, one can conclude that the existence of non-radiating 
spinning solitons, although not forbidden, is not guaranteed either. Even for known static 
solitons their spinning generalizations may or may not exist, and if they do exist, this should 
be considered as something exceptional rather than the general rule. At first view, such a 
conclusion may seem to contradict our experience, since normally one knows that compact 
objects can spin. However, we are not saying that generic field theory solitons cannot spin. 
They always can, but it seems that they should generically radiate at the same time. One can 
always 'give a kick' to any static soliton, as for example to the magnetic monopole or to any 
of the knotted solitons of Faddeev-Skyrme shown in Figj6l so that they will start spinning. 
However, at the same time they will start radiating away all the received energy and angular 
momentum, till they relax back to the original static, non-spinning configuration. It seems 
that 'spinning and radiating' represents the generic way the field systems behave, while 
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'spinning without radiating' should rather be considered as something exceptional, possible 
only in some special field theory models. 

Perhaps the best way to really establish the existence of spinning and non-radiating 
solitons would be to construct them as solutions of the elliptic boundary value problem 
obtained by separating the time variable. In gauge field theory, as for example in the local 



U(l)xU(l) Witten's model |l77j ]. one could look for manifestly stationary and axisymmetric 
vortons, in which case there is no radiation. However, since the angular momentum in 
this case can be expressed by the surface integral fl4.66p and so is determined only by the 
asymptotic behavior of the fields, there are higher chances that it could vanish. For global 
field theories one could consider systems with non-manifest symmetries, in which cases there 
are better chances to have a non-zero angular momentum, but also higher chances to have 
radiation, unless the no-radiation condition fl4.7ip is fulfilled. 

V. EXPLICIT EXAMPLES OF STATIONARY SPINNING SOLITONS 

Before coming to vortons, one can wonder, in view of the above discussion, if there exist 
at all any known examples of non-radiating spinning solitons. As we said, for the SU(2) 
magnetic monopoles, dyons and sphalerons there are no spinning generalizations, at least in 
the manifestly stationary and manifestly axisymmetric case. Nevertheless, explicit examples 
of spinning solitons in Minkowski space in 3 + 1 dimensions exist and below we shall review 
all known solutions of this type. These are the spinning Q-balls, spinning Skyrmions and 
also rotating monopole-antimonopole pairs, apart from the vortons. Interestingly, it seems 
that spinning solitons are in some sense more easily constructed in curved space, because the 
spinning degrees of freedom can be naturally associated to the non-radiative dipole moment 
of gravitational field. For this reason a number of articles cited below actually describe the 
Minkowski space spinning solitons only as a special limit of the more general, self-gravitating 
configurations. 

A. g-balls 

This is an important for our discussion example which shares many features with the 
vortons. At the same time, spinning Q-balls are simpler than vortons, and so they can be 



42 

used to introduce a number of notions to be applied later. We shall therefore discuss them 
in some detail. 



Q-balls have been introduced by Coleman 39|]. These are non-topological solitons 116 1 



found in a theory with a single complex scalar field with the Lagrangian density 

/:Q[$]=9^$*9^$-f/(|<D|), (5.72) 

the corresponding field equation being 

^m5^^ + ^^ = 0, (5.73) 

while the energy-momentum tensor 

T^u = 9^$*9,$ + 9.$*9^$ - q^Xq- (5.74) 

The potential U should have the absolute minimum, ^(0) = 0, and should also satisfy the 

condition 

TT(f\ _ 1 rl'^TT 

(5.75) 



uj = mm — — — < LO — 



/=o 



7" P + 2 dp 

whose meaning will be explained below. The value of u^ determines the mass of the field 
quanta, M = a;+. For the condition fl5.75p to be fulfilled U should (if it is even in |$|) contain 
powers of |<l>| higher than four, which means that the theory cannot be renormalizable 39 1. 
A convenient choice is 

[/ = A|$p(|$|^-a|<l>|2 + 6), (5.76) 

where A, a, b are positive constants, so that c^i = Xb and a;^ = u;i(l — a'^/Ab). In our 



numerics below we shall always choose A = l,a = 2,6 = l.l 168l ]. which is not a restriction 
(for u^ 7^ 0) since A, a can be changed by rescaling the coordinates and field, while the 
mass M = y/Xb enters the field equation (15.831) only in the combination M^ — cu'^ where cu is 
another free parameter. 

The global invariance of the theory under $ -^ $e*" implies the conservation of the 
Noether charge 

Q = I Udt^*^ - ^*dt^)d^x . (5.77) 

The scaling argument of Derrick applies for the theory (I5.72p and rules out all static solutions 
with finite energy. Therefore, in order to circumvent this argument, solutions should depend 
on time, 

$ = 0(x)e^'^*, (5.78) 
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in which case the Noether charge is 

Q = 2u (p^d^x = 2ujJ\f. (5.79) 

In what follows we shall assume u to be positive. More explicitly, with 

E2= j{V(t)fd^^, Eo=fud^x, (5.80) 

the Lagrangian reads 



CQd^^ = uj^M -Eq-E2. (5.81) 

Under scale transformations x -^ Ax one has L — > K^{uj'^J\f — Eq) — AE2 so that L will be 
stationary for A = 1 if the virial relation is fulfilled, 

3uJ^^f = 3Eo + E2, (5.82) 

which is only possible (for Eq 7^ 0, £'2 7^ 0) if cj 7^ 0. 

Q-balls are finite energy solutions of the Lagrangian field equations for the ansatz (15.781) . 

{A+uj^)(j)=^<p. (5.83) 



Even though $ depends on time, T^j/ is time independent - the system is non-manifestly 
stationary. Equivalently, Q-balls can be obtained by minimizing the total energy 

E= f T^d^x = /(c^VP + I V0P + U)d^x = uo^U + E0 + E2 (5.84) 

by keeping fixed the charge Q. Indeed, rewriting the energy as 

E=^ + Eo + E2 (5.85) 

shows that minimizing E with Q fixed is equivalent to extremizing L, since 

5E = -^ 5^f + 5{Eo + E2) = -ujHM + 5iEo + E2) = -5L , (5.86) 

so that the on-shell condition 6L = follows from SE = 0. It is also instructive to see how 
the same thing comes about within the Lagrange multiplier method. Introducing 

,2|i|2 , iv7J,|2 , rrNj3„ I ,. /o, , / Ij,|2j3 



Eq = / {lu'\(I)\' + \V(I)\' + f/)d'x + i2[2iu 101 Vx - Q , (5.87) 
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the condition dEq/dn = fixes the charge, the condition dEg/du = gives /i = —u 
insuring that Eq = L + const., so that the condition 6Eq/6(J) = reproduces the field 
equation. 

Equivalently, one can minimize the potential energy of the field, Eq + E2, by keeping fixed 



^f= l^pd^x (5.88) 

via extremizing the functional 

E^ = /"(I V0I' + f/)rf'x + fii( f 101 Vx - AT j . (5.89) 

The condition dEj^/dfii = imposes the constraint (I5.88p . while the condition 5Ej^/5(j) = 
reproduces the field equation (I5.83P with u"^ = —fii. The time dependence of the solutions 
is very implicit in this approach, since u appears only as a Lagrange multiplier. It follows 

from Eq. (15:861) that 

_ djEo + E,) 
^ - g^ , (5.90) 

so that fixing A/" fixes also the charge Q = 2uj{N')M . In the non-relativistic theory M can 
be viewed as the particle number (see Sec JVII A] below) . 

Let us also note that using the charge definition (I5.79P and the virial relation (I5.82p the 
energy 05.841) can be rewritten as 

E = uQ + -E2. (5.91) 

o 



Q-balls could exist in the supersymrnetric extensions of Standard Model 113l | and could 



perhaps contribute to the dark matter 114 1 



1. Non-spinning Q-balls 

since this will help 



Let us briefiy consider the simplest spherically symmetric Q-balls ^ 
to understand more complex, spinning solutions. Setting $ = e*'^/(r), the real amplitude 
/(r) satisfies the equation 

r + -/' + -V = ^^. (5.92) 

r 2 oj 

For the energy to be finite / should vanish at infinity, so that for large r one has 

/ ~ - exp{-VM2-cu2r}. (5.93) 
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FIG. 11: Left: the amplitude /(r) for the n = 0, 1,2 spherically symmetric Q-bahs with uj = 0.2. 
Right: the effective potential C/efr in Eq. (j5.94p for three values of uj. 



Solutions of Eg. (15. 921) c omp rise an infinite family labeled by an integer n 



0,1,... 



counting the nodes of /(r) 168| (see Fig JTTll. The energy increases with n. These solutions 



admit a simple qualitative interpretation 



39| . since Eq. (l5.92p is then equivalent to 
dr 



r + U,M=E-A / / 



f/2 



(5.94) 



This describes a particle moving with friction in one-dimensional potential Ucs{f) = ^"^P — 
U{f) shown in FigITT| the integration constant E playing a role of the total energy. At 
the 'moment' r = the 'particle' rests at a point with coordinate /o close to the potential 
maximum A (see Fig JTTI) . so that its energy is E = Ucs{fo)- Then it starts moving to the 
left, dissipating its energy as it goes. One can adjust the value of /o so that for r — > oo it 
dissipates all its energy and arrives at the local maximum B of the potential with zero total 
energy to rest there - either directly or after n oscillations between the two potential hills. 

It follows that UJ should belong to the interval u;_ < u < u^, the condition Eq. (l5.75p 
making sure that this interval is non-empty. Q-balls become large as cj — > u;±, their charge 
and energy growing without bounds. 

As UJ -^ UJ- the maximum A of the potential descends towards the position D. The 
'particle' then stays for a 'long time' at D, till the friction term becomes suppressed by 
the 1/r factor. Then it starts moving, crosses the potential well in a finite 'time' Ar and 
asymptotically approaches the potential maximum B. In the simplest n = case solutions 
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in this limit can be described by a smoothed step-function, 

f{r)^foQ{R-r), (5.95) 

where R ^ oo as uj -^ u^ and /' ^ /o/Ar in a region of a fixed size Ar ^ i? around r = R. 
This is sometimes called Q-balls in the thin wall approximation 39|. Their charge Q ^ R^ 



while E2 ~ R^ so that for large R one can neglect the second term in Eq. (l5.9ip . which 
gives for their energy E = uj^Q. The same value is obtained inserting (15.951) to Eg. (15.85 1 ) . 
neglecting the E2 term and minimizing with respect to -R, /q. 

When uo — > uo^ the local minimum C of the potential approaches the local maximum B. 
The potential in their vicinity becomes approximately Ucs ~ f^{f'^—M^) with M^, = lJ\—lJ^. 



Inserting this to Eq. (l5.94p gives Q-balls in the thick wall approximation 112 | 



f{r)^^PlM^y{M^r), (5.96) 

where y{x) fulfills y" + {2/x)y' + (y^ — l)y = 0. Solutions of this equation satisfy y{0) = yn 
and y ~ e~^ for large x (one has yo = 4.33). Although / ~ M^^ is small in this case, the 
typical configuration radius R ~ 1/M^ is large. The charge is also large, Q ~ 1/M^, while 
E2 ~ M(^, so that one can again neglect the second term in Eq. (l5.9ip . which gives for the 
energy E = uj^Q. 

2. Spinning Q-halls 
The spinning, axially symmetric generaliz ation s for the spherically symmetric Q-bal 



lave b een c onstructed by Volkov and Wohnert 168l | and have been analysed also in Refs. [10|] , 



102| . |l03l ]. [301 ]. These non- manifestly stationary and non-manifestly axisymmetric solu- 
tions provide the first explicitly constructed example of stationary, spinning solitons in 
Minkowski space in a relativistic field theory in 3 + 1 spacetime dimensions. They have the 
typical ring structure, very similar to that for vortons, and they can also have the same 
topology as vortons. They admit simple generalizations to the gauged case. Spinning Q- 
balls thus provide a simple prototype example of vortons. Their analogs also exist in the 



context of non-linear optics 1251 ]. where they describe spinning light pulses (see Sec JVII B] 
below) . 
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Spinning Q-balls have a spinning phase: the field being non-manifestly stationary and 
non-manifestly axisymmetric, 

q, = ^i{^t+mv)f^r,^). (5.97) 

The real field amphtude f{r,'d) satisfies the equation 

For m ^ the energy will be finite if only / vanishes at the symmetry axis, while at infinity 
the asymptotic behavior f l5.93p still applies. The energy momentum T^^ depends only on r, -d 
and the angular momentum is 

J = T^d^x = 2muj / /^t/^x = 2mcj7V, (5.99) 

so that it is classically quantized as 

J = mQ. (5.100) 

In view of this relation, spinning Q-balls correspond to minima of energy with fixed angular 
momentum. They could therefore be obtained by extremizing 

with constant J,m, the corresponding extremum condition being given by Eq. (15.981) . A 
direct minimization of this functional carried out in Ref. [lO| (although for a different choice 
of the potential U) suggests that its extrema exist but seem to have negative directions. 
One should therefore integrate the field equation to construct these solutions. 

Taking into account the boundary conditions at the 2;-axis and at infinity, Eq. (I5.98P with 
m 7^ admi ts two different types of solutions determined by the behaviour of / under 



z^ -z [l68|. 



For even-parity solutions, called tti"*", the amplitude /(r, '&), energy-momentum and charge 
densities are maximal in the equatorial plane and the energy is concentrated in a toroidal 
region encircling the 2;-axis (see Fig JT2l) . 

For odd-parity solutions, called m~ , the amplitude f{r^d) vanishes in the equatorial 
plane, while the energy-momentum and charge densities show two maxima located symmet- 
rically with respect to the plane, so that the solutions exhibit in this case a double torus or 
dumbbell-like structure (see Fig lTSl) . 
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FIG. 12: The amplitude /(r, ??), energy density Tq and the energy isosurface with Tq = 0.65 for 
the 1+ spinning Q-ball with uj = 0.9. 







FIG. 13: The amplitude /(r, ??), energy density Tq and the energy isosurface with Tq = 0.65 for 
the 1~ spinning Q-ball with to = 0.9. 



Apart from choosing the parity and the value of m, solving the differential equation fl5.98p 
also requires choosing u^ as an input parameter, the energy and charge being computed from 
the numerical output. As in the spherically symmetric case, spinning solutions are also found 
to exist only in a finite frequency range, uj^ < u < u;+. Here U+ = M as in Eq.f l5.75l) . while 
U- seems to be ?7i-dependent. As u approaches the limiting values, the energy and charge 
seem to grow without bounds (see Fig lT^ . while somewhere in between there is a critical 
value of frequency, cUcrit, for whic h both E{uj) and Q{uj) attain their minimal values. This 
behaviour is found for all m |l02l | and it is qualitatively the same as in the m = case. 



49 

The existence of a minimal value of the Q-ball charge implies that the angular momentum 
(15.1001) cannot be arbitrarily small. Q-balls cannot therefore rotate slowly. They show a 
discrete spectrum of spinning excitations. 

Although it is difficult to qualitatively analyse the behaviour of the m ^ solutions, 
there are some analogies with the m = case. Specifically, Q-balls become large as u; — ^ uj±. 
For uj -^ uj^ they can be viewed as squashed spheroids, homogeneously filled inside, which 
reminds of the thin wall approximation. Unfortunately, an analog of the step function 
solution fl5.95p does not directly apply for m ^ 0, since the function / cannot be constant 
inside the ball because it must vanish at the symmetry axis. Instead, / increases as one moves 
away from the axis and reaches maximal values at the surface of the spheroids, after which 
it rapidly goes to zero. The energy density is approximately constant inside the spheroid, 
with a slight increase at its surface, and rapidly vanishes outside it. For u; — > u;+ solutions 
also become large spheroids, but this time they are hollow, with the maximal energy density 
concentrated at the surface and being close to zero everywhere else. It seems that the thick 
wall approximation (15.961) can be directly generalized to the axially symmetric case to give 
/(r, 6) ^ V2M^ y{M^r, 6) where M^ ^ and y{r, 6) fulfills (A - ^^^ + y^ - l)y = 0. 

If one considers Q and not u as the solution parameter then, using E{uj) and Q{uj) to 
express E in terms of Q , one discovers that the function E{Q) is double-valued with a cusp. 



as shown in Fig{Tl] 102| . For a given Q there are thus two different spinning Q-ball solutions 
with different energies. Solutions from the less energetic branch correspond to the uj < cUcrit 
parts of the E^u), Q{uj) curves in FigHH As in the m = case, the plots shown in FigJT^ 
demonstrate for large Q the linear dependence, E = uj±Q. This can again be explained by 
the general relation (I5.9ip . implying that for large Q one can neglect the gradient energy E2 
as compared to Q. 

The existence of two different solutions with the same Q suggests that the more energetic 
of them is unstable. In fact, it follows from Eq. (l5.9ip that the energy-to- charge ratio for 
Q-balls is E IQ = u; + 2E2/3Q. For the upper branch solutions one has u -^ uj+ = M for 
large Q, implying that E/Q = M + 0{Q~^) where the subleading term can in principle 
be positive or negative, depending on the details of the Q-ball potential. In the present 
case it is positive, as can be seen in FigJT4l so that the upper branch Q-balls are unstable 
with respect to decay into free particles. The same argument for the lower branch solutions 
gives E/Q = uo^ < M for large Q, so that they cannot decay into free particles. However, 
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FIG. 14: Energy and charge E{uj),Q{uj) (left) and E{Q) (right) for the 1^ spinning Q-balls. For 
large Q one has E{Q) = uj^Q and E[Q) = uj^Q for the lower an upper branch, respectively. The 
upper branch solutions are unstable, since E/{MQ) > 1, while the lower branch solutions are 
stable for large enough Q. 



there can be other decay modes, and so the stability analysis is needed. In fact, the m = 
lower branch Q-balls are known to be stable 4^. The stability analysis of the even parity 
solutions with m ^ will be described below in Sec JVIIBl It seems that for |?7i| = 1 the 
lower branch solutions are stable only for large enough Q, while already for |??7,| = 2 all of 
them are unstable - they seem to decay by splitting into several non-spinning Q-balls. 

It is plausible that spinning generalizations could also be constructed for the excited 
spherically symmetric Q-balls with n > 0. The complete family of spinning Q-balls should 
therefore contain not only m^ solutions, but also the excited {n,m^) solutions with n = 
1,2, .. . for which the amplitude f{r,§) exhibits nodes. More precisely, parametrizing the 
(p, z) plane by a complex variable w = p + iz, one can expect the amplitude f{w) of 
the {n,m~^) and {n,m~) solutions to have the same zeros as, respectively, ^{F^{w)) and 
Q{F~{w)), where 

n 

F^{w) = llaf{w-pf), (5.102) 

i=i 

with some suitably chosen af E C and pf > 0. However, very little is known about such 
excited solutions at present. 
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3. Twisted Q -halls 

Here and in the next subsection we briefly describe our new results on further general- 
izations of the spinning Q-balls, not yet discussed in the literature. Let us consider again 
the theory fl5.72p . but generalize the fleld ansatz to include an independent phase, 



$ _ pi{(^t+m<p-nip{r,9)) 



fir, 9) = {X{r,^)+iY{r,^))e 



i{iL!t+mip) 



(5.103) 



We require the phase function il){r,'d) to increase by 27r after one revolution around the 
contour C shown in Fig|2l The overall phase ut + mip — nilj{r,6) then winds around the 
circle 5* and along the contour C, exactly as for the Faddeev-Hopf field Eg. (12. 101) . For regular 
fields /(r, ■&) vanishes at C and has n zeros inside C. The solutions are thus characterized by 
two integers (n, m) giving rise to the 'topological charge' A^ = nm, although now this does 
not represent a genuine topological invariant. We shall call Q-balls with A^ 7^ twisted, 
by analogy with twisted loops in the Faddeev-Hopf theory, while those with A^ = n = 
described above will be called simply spinning or 'non- twisted'. 





FIG. 15: Profiles of the (1, 1) twisted Q-ball with uj = 0.9. 



The field equations read 

A- 
A- 



m 



m 



^ 2 +uj'^]X 
r^ sin {} j 2 



1 dU{^PO + Y^) 



dX 



uj^\Y 



1 dU{^X^ + F2) 



(5.104) 



r2 sin^ ^ ' " j ' 2 dY 

where one can require X, Y to be symmetric and antisymmetric, respectively, with respect to 
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reflections in the equatorial plane. The simplest twisted solutions are obtained for n = m = 1 
and are shown in FiglTSl 





FIG. 16: Energy density for the (1, 1) twisted Q-ball with uj = 0.9 (left) and the isosurfaces of 
constant energy density with Tq = 0.3 (right). 

The proflles of the energy density look spectacular; see FigJT6l It is unclear at present 
whether these solutions are dynamically stable, but they certainly exist as solutions of the 
elliptic system (I5.104p . The twisted Q-balls are much more heavy than the non- twisted ones. 
For example, for uj = 0.9 we found for the twisted (1, 1) solution the energy E = 929.34, while 
the non-twisted 1+ and 1~ spinning Q-bahs have, respectively, E = 169.27 and E = 338.66. 

4- Spinning gauged Q -balls 



Another possibility to generalize the spinning Q-balls is to couple them to a gauge fleld 
within the theory 

£gQ[A^,$] = -1f^,F^^ + £q[$]. (5.105) 

Here F^^, = dfj_Ai, — d^A^ while £q[$] is the same Q-ball Lagrangian as in (I5.72p but with 
the derivatives of $ replaced by the covariant derivatives, (9^$ -^ -D^$ = (9^ — igA^j)^ 
where g is the gauge coupling constant. This theory is in fact the U(l) version of the general 
gauge fleld model (I3.50p . The local gauge transformations (13.521) now read $ -^ $e*^". 
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Afj_ ^ A^ + d^a while the field equations (13.531) assume the form 

D,D^$ = -^$. (5.106) 

The conserved Noether charge analogues to the Q-ball charge Eg. (15.771) is 

Q= fj^d'^=-<isdS = ^^, (5.107) 

where Qei is the electric charge. Spherically symmetric solutions of this model were discussed 



in Ref. 



115 |. 



We make the ansatz 

Af,dx^' = Ao{r,^)dt + A^{r, ^) sin ed^, $ = /(r, ^)e*('"^+'^*), (5.108) 

with real /(r, ■(9), and require the gauge field to vanish at infinity, while at the 2- axis A^ = 
dgAQ = 0. Although $ depends on t, ip, this dependence can be gauged away, so that the 
system is manifestly stationary and manifestly axially symmetric. Numerically solving the 
field equations gives spinning Q-balls with a long range gauge field which behaves for large 
r as (m being the magnetic dipole moment) 

Qci . msin6' , _,^^, 

^0 = ^^ + • • • , A^ = ^— + . . . . (5.109) 

In the limit (7 — > they reduce to the non-gauged spinning Q-balls. It seems that solutions 
exist if only g does not exceed a certain maximal value gmaxi^) (see FigJTTT). This feature 
can be understood qualitatively 1151 ]: since Q-balls can be viewed as condensate states of 
mutually attracting scalar particles, gauging them creates an electric repulsion that destroys 
the condensate for large enough g. 

Since the symmetries of the solutions are manifest, one can use the surface integral 
formula (14.661) for the angular momentum. Applying the symmetry equations (I4.64p to the 
ansatz (15.1081) gives W^ = m/g, inserting which to (14.661) and using (I5.109P yields 

■J = f{j-\) £dS = ^^ = Qm. (5.110) 

This shows that spinning is possible in manifestly stationary and manifestly axially sym- 
metric systems where all the spinning phases can be gauged away. 
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FIG. 17: Left: the energy E{uj) and E[Q)/[QM) for several values oi g for the l"*" gauged Q-balls. 
Right: the energy E{g) for several values of lo for the 1"^ solutions. The plots of Q{uj) and Q{g) 
look qualitatively similar. 

As seen in FigJTTt the dependence of solutions on uj is similar to that in the ungauged 
case: they exist for a limited range of u. Expressing E in terms of Q gives again a two- 
branch function E{Q), the solutions from the lower branch being stable for large enough Q 
with respect to decay into free particles. E, Q cannot be arbitrarily large for g ^ 0, since the 
electric and scalar field contributions to the energy grow as Ql^ ~ Q^ and Q, respectively, 
and so for large Q the electric term dominates and destroys the soliton. 

For small g solutions can be represented as $ = $'-°'' + g"^^^"^^ + . . . and A^ = gAii + . . . 
where ^^^^ is the non-gauged Q-ball. The energy is E{g) = E^^^ + g'^ E^'^^ + . . . and calculating 
^(2) reveals that it is not sign definite - due to the electric field contribution. The plots in 
FigHTI show that it can be both positive and negative, depending on the solution. 

5. Spinning interacting Q-balls 



Yet another way to generalize Q-balls is to couple to each other several copies of the the- 
ory. This allows one to consider 'non-linear superpositions' of individual Q-balls. Although 
this idea has already been considered in the literature 30|], we have reconsidered it and 
found some curious solutions which could be interesting in the context of our discussion. In 
the simplest case one can choose 



L[$i, $2] = lq{<^i) + lq($2) - 7i*iri*2r' 



(5.111) 
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FIG. 18: Energy isosurfaces with Tq = 0.1 for the mi = 0, m2 = 3 'Saturn' (left) and with 
Tq = 0.3 for the mi = 1, m2 = 3 bi-ring (right) solutions in the theory (jS.llip . 

where Lq(<I>i) and Lq($2) are two copies of the Q-ball Lagrangian (I5.72p . Setting 

$^ _ ^i{c.it+m,^) j^^^^ Q^^ ^^ ^ ^i{^2t+m,^) f^^^^ Q^^ (5.112) 

it is interesting to consider solutions with rrii ^ 1712- Choosing ct^i = u;2 = 1, 7 = 1 and 




▲ 



K 




FIG. 19: The energy density Tq and the Tq = 0.1 energy isosurface for the ttii = 0, m2 = 3 'hoop' 
solution in the theory ()5.11ip . 

restricting to the even parity sector, we find solutions of the 'Saturn type', with a central 
concentration of the energy density produced by the first scalar with mi = and surrounded 
by a ring created by the second scalar with m2 > (Fig JTSjl . solutions for 7712 > mi > with 
bi-ring profiles (Fig JTSj) . as well as solutions describing a superposition of a straight 'Q-ball 
vortex' with mi = 'hooped' by a Q-ball ring with m,2 > (Fig JTOl) . 

B. Skyrmions 

The first known example of topological solitons in 3 + 1 dimensions was suggested almost 
50 years ago within the non-linear relativistic field theory model proposed by T.H.R. Skyrme 
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157|. These solitons are now called skyrmions. Skyrme himself considered them as field 
theoretic realizations of baryons. Nowadays the Skyrme model is regarded as an effective, 
low energy approximation of QCD |5|,[j]. In this approximation the static skyrmions are 
promoted to spinning objects by making use of the effective rigid body approximation that 
will be considered below. The masses of the spinning skyrmions obtained in this way are 
then compared to the hadron masses. 

The question of whether spinning skyrmions really exist as stationary field theory objects 
was addressed only very recently in Ref. [l8|, whose authors constructed spinning skyrmions 
by applying the same mechanism as for the Q-balls and arrived at conclusions which differ 
considerably from those obtained within the rigid body approximation. 

Since there exist excellent descriptions of the Skyrme model and its solutions in the 



literature 120| , 12 1| , we shall very briefly summarize only the features essential for our 
discussion. The fundamental fleld variables in the theory comprise an SU(2)-valued matrix 
U{x^) which satisfles fleld equations for the Lagrangian 

Cs[U] = tr r^d^U^d>'U+^ [d,U\dMWU^,d'U]\ . (5.113) 

It is convenient to introduce a pure gauge SU(2) connection, 

A,, = lU^d^U = r'^Al, d^A, - d,A^ = i[A^, A,], (5.114) 

in terms of which the fleld equations read 

d'^ (a^ - i [[A^, A.], A'']] = 0. (5.115) 

The energy for static flelds is 

E[U] = jtr Q(A,)^ - ^([A„ A,])^) d^^. (5.116) 

For the energy to be flnite the fleld f/(x) should approach a constant value at inflnity, which 
can be chosen to be the unit matrix, so that lim|x|--+oo U{^) = 1. This allows one to replace 
M^ by its one-point compactiflcation S^. Since SU(2) is topologically also S^ it follows that 
any flnite energy fleld conflguration can be viewed as a map, f/(x) : S^ -^ S^, and can 
therefore be characterized by the integer degree of map, 

N[U] = ^ /" tr (e.jkAiAjAk) d'x . (5.117) 
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This topological charge is called in Skyrme theory 'baryon number'. Fixing it the energy 
obeys the Bogomornyi type inequality 

E[U] > 12V2n'^\N[U]\, (5.118) 

which can be easily obtained by rearranging terms in the integrand in f l5.116p . The existence 
of this lower bound suggests looking for energy minimizers in each topological sector. 

It will be convenient for what follows to parametrize the matrix U in terms of two complex 
scalar fields 0, a subject to the normalization condition |0p + |crp = 1 as 

f/= r" . (5.119) 

yta <p* ) 

Simplest skyrmions are spherically symmetric, U = exp(ix(r) x ■ T/r), hence 

= cosx(r) + 2sinx(r) COS"!?, a = sinx(r) sin'i^e*'^. (5.120) 

Inserting this to Eqs. (l5.115l) the variables separate and the equations reduce to an ODE 
for x(^)- This equation admits globally regular solutions with the boundary conditions 
x(0) = A^vr and x(oo) = for any value of the topological charge A^. However, it seems that 
only for A^ = ±1 these solutions correspond to the absolute energy minimum, while those 
for |A^| > 1 are local minima or saddle points. Global energy minima for |A^| > 1 are not 
spherically symmetric and can be obtained by directly minimizing the energy (I5.116p . For 
|A^| = 2 they are axially symmetric, with the constant baryon number isosurfaces having 
toroidal shape, so they are somewhat similar to knots. However, higher charge skyrmions 
do not resemble knots at all and look like polyhedral shells. Their detailed description can 



be found in the monograph 121 |. 



1. Skyrme versus Faddeev- Skyrme 

It is worth emphasizing the similarity between the Skyrme and Faddeev- Skyrme models. 
Both can be represented in the form (up to normalization) 

£[0] = d,cl> ■ 9^0 - ^ (d.cp ■ d'^cj.f + ^ {d^cj, ■ 9,0)(9^0 ■ d'^cj,) (5.121) 

where (^ is a vector normalized by the condition </> ■ = 1. In the Faddeev-Skyrme case 
one has cp = (p" with a = 1,2,3 so that the target space is S"^. In the Skyrme model 



58 

one has (f> = (0°, 0") with a = 1, 2, 3 so that the target space is 5*^, the Skyrme field being 
U = (f)^+i(f)°'T°-. One can also construct a generalized theory whose t arget space 'interpolates' 



between 5*^ and S^ and so it includes both models as special cases [175I ]. 

The Faddeev-Skyrme model can be regarded as a consistent truncation of the Skyrme 
theory. Specifically, if n = ra" is a solution of the Faddeev-Skyrme model, then = (0, n") 
solves the equation s of the Skyrme model. As a result, knot solitons can be embedded into 



the Skyrme model 123|, 37|], although their stability properties will then be different. 

Another embedding of knot solitons into Skyrme theory is obtained by expressing the 
fields Afj, and n in terms of the CP^ coordinates (f>,a using Eqs.f l2.28l) .f l2.30l) .f l2.3ip . The 
Skyrme field U is then given by Eq. (15. 1191) . Although this correspondence does not map 
solutions to solutions, the advantage now is that the baryon number N[U] is exactly equal 
to the Hopf charge A^[n]. This follows from the fact that Aj^ = —A^ and jF^j^ = 2(A)jA^ — 
A;',A^) due to Eq. (15.1140 . With this, the expression (15.1171) for the baryon number exactly 
reduces to the expression ( 12. 9p for the Hopf charge. For example, using the parametrization 
(I2.32P for the CP^ scalars 0, a with the Hopf charge A^[n] = nm gives with (I5.119P the 
Skyrme field U with N[U] = nm. 

The Hopf charge A^[n] will determine in this case also the winding number of the pure 
gauge fi eld A ^. = iWdkU. If n is an energy minimizer then the pure gauge A^ is maximally 



abelian 164 1 



The correspondence can also be used in the opposite direction to construct n with a 
given Hopf charge out of U with a given baryon number. If U is spherically symmetric, 
then n is axially symmetric, as can be seen comparing Eq. (l5.120p and Eq. (l2.16p . Using 



approximations for U gives in this way approximate solutions for n 173 1. 



2. Spinning skyrmions 

Let us consider the A^ = 1 static, spherically symmetric skyrmion described by Eq. (l5.120p . 
Its spinning analog in the rigid body approximation is obtained by simply replacing in 
Eq.f l5:T20|) 

if^if + ut, (5.122) 

which gives a non-zero value to the angular momentum. It is in fact precisely this type of 
approximation that is often implicitly assumed in the literature when talking about spinning 
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solitons. In this approximation spinning solitons do not change their shape and the radiation 
effects are neglected, which are natural assumptions if the rotation is slow. However, if uj is 
not small, then this description does not approximate the true spinning solutions any more 
but simply gives field configurations with J 7^ 0. If one uses them as initial data, then their 
temporal evolutions will be certainly accompanied by radiation carrying J away, and it is 
not clear if the system will finally relax to a state with J 7^ 0. 

In order to construct truly spinning skyrmions with A^ = 1, Battye, Krusch and Sut- 



cliffe (BKS) 18| generalize the spherically symmetric ansatz fl5.119p . fl5.120p to the axially 



symmetric, non-manifestly stationary one, which can be parametrized as 

= cos^e^'^, a = sin^e*(^+'^*). (5.123) 

It is instructive to compare this with the ansatz fl2.32p for the axially symmetric hopfions 
with m = n = 1. They also add to the Lagrangian f l5.113p the mass term 

M^tr(f/-1). (5.124) 

Instead of solving the field equations for 0(r, '&), ip{r, '&) they minimize the energy with fixed 
J, whose expression is similar to the one in Eq. fl5.101|) . The frequency u for their solutions 
is a parameter in the interval 

0<u^ < M^, (5.125) 

both the energy and angular momentum increasing with u. The latter feature is very 
interesting - for spinning skyrmions J is a continuous parameter that can be arbitrarily 
small, so that they can rotate slowly. To compare, Q-balls cannot rotate slowly. Although 
BKS do not emphasize this, it seems that both E and J should blow up as u"^ -^ M^ since 



the fields localized by the factor exp(— vAP^^u^r) become then long range. 

Trying to adjust M^ and the energy scale to reproduce the pion and hadron masses and 
their spins, BKS discover that this requires moving to the parameter region where deviations 
of spinning skyrmions from spherical symmetry are large. They conclude that the rigid body 
approximation does not provide an adequate description of spinning skyrmions. 

The results of BKS have been confirmed by solving the field equations by loannidou. 



Kleihaus and Kunz 



94| (who studied in fact gravity-coupled skyrmions, but considered the 



Minkowski space limit as well). 
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The spinning skyrmions have been obtained by BKS by minimizing the energy within the 
axially symmetric ansatz (15.1231) . Strictly speaking, this does not exclude the possibility of 
non-axially symmetric instabilities. However, since the non-spinning skyrmions are stable, 
one can expect that at least for small J their spinning analogs should be stable as well. 

Spinning solitons have also been studied in the context of the baby Skyrme model 



133| .|25l]. This is simply the Faddeev- Skyrme model restricted to 2 + 1 spacetime dimensions 
and written in the form fl5.12ip with an additional mass term M^(0^ ~ !)• Spinning fields 
are then chosen to be 

0i + Z02 _ sine(p)e^('^*+'^), cf)^ = cose(p). (5.126) 

The field equations reduce in this case to an ODE for 6(p) whose solutions are easy to study. 
Solutions exist for cu < M and both E and J blow up as u; — ^ M. 

Trying to further increase u, solutions become oscillatory, since one has O ~ 



exp(— vAP^^^p) for large p, both E and J being then infinite. However, one can consider 
the initial data of the form ( ]5.126p . where u > M and 0(p) vanishes identically for large p. 
E, J will then be finite. Evolving these data in time, the system radiates away a fraction of 
its ener gy an d angular momentum and relaxes to a stationary, rotating configuration with 
Lj < M (l33|. 



3. Spinning gauged skyrmions 
Spinning solutions have also been constructed within the gauged version of the Skyrme 



model 3l|] by Radu and Tchrakian 138 



These solutions are manifestly stationary and 
manifestly axisymmetric. Parametrizing the Skyrme field U in terms of two complex scalars 
(f), a according to Eg. (15. 1191) . the gauged Skyrme model is obtained from the Skyrme model 
as 



C,s[A^,4>,a] = -- E.^E'^'' + Cs[^,a]. (5.127) 



1 
4 
Here Cs is the Skyrme Lagrangian (15.1131) with the mass term (15.1241) included, E^^, = 

d^A^ — di^A^, and the derivatives of the field a in Cs[4>.,o'] are replaced by the covariant 

derivatives, d^a -^ {d^ — igA^)a, while the derivatives (p do not change, so that the theory 

is invariant under local U(l) gauge transformations 

A^^A^ + d^a, a^e^'a, 0^0. (5.128) 
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Radu and Tchrakian 



138l ] generalize the ansatz (15.1231) as 



A^,dx^' = AQdt + A^suidd^, (/. = cos^e^'^, a = sin ^ e'("'^+^*) , (5.129) 

where Ao^A^p^Q^ip depend on r^d. This ansatz is manifestly stationary and manifestly 
axially symmetric, since its t, (/9-dependence can be removed by a gauge transformation. 
The topological charge fl5.117p is A^ = m. Integrating the field equations gives globally 
regular solutions with finite energy and a longrange gauge field, which is at large r 

Ao = — + ..., A^ = ^ sin^ + ..., (5.130) 

where Qei, m are the electric charge and dipole moment. This feature is similar to that for 
the gauged Q-balls in Eg. (15. 1091) . The angular momentum can be calculated in the same 
way as in Eg. (15. 1101) . with a similar result, 

j=^Jl9^. (5.131) 

However, unlike for Q-balls, J can be arbitrarily small if uj is small. The dependence of 
solutions on uo is gualitatively similar to that for the global skyrmions: they exist only for 
a finite freguency range, both E and J growing with u. For cj = the solutions reduce to 



the non-spinning gauged skyrmions found in Ref. [1341 ]. 

The stability of spinning gauged skyrmions has not been studied. They probably have 
less chance to be stable than their global counterparts - since gauging introduces additional 
degrees of freedom that can produce instabilities. 

C. Rotating monopole-antimonopole pairs 

There are other known solutions with non-zero angular momentum. They are mani- 
festly stationary and manifestly axisymmetric, but they describe rotations in multisoliton 
systems and not spinning of a single soliton. Let us consider the already mentioned above 
monopole-antimonopole s oluti ons in the Yang-Mills-Higgs theory (13.501) . Their existence was 



demonstrated by Taubes [l61l | and they were explicitly constructed by Kleihaus and Kunz 
100| within the ansatz (I3.55P with k = 2, m = 1. These solutions are static, purely mag- 
netic and manifestly axisymmetric (since their (/^-dependence can be gauged away) and with 
vanishing magnetic charge. The Higgs field $ shows two simple zeros at two spatial points 
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separated by a finite distance (see Fig j20|) . These two points correspond to the positions of 
the monopole and antimonopole. The analysis of the charge and current distributions similar 
to that described in Sec JIIIBl shows that the system also contains a circular electric current, 
as shown in Figl20l The magnetic field created by this current acts against the Coulom- 
bian attraction between the monopole and antimonopole, which presumably stabilizes the 
system. 





FIG. 20: Left: the Higgs field amplitude for the nionopole-antinionopole solution in the limit where 
the Higgs field potential vanishes. It is shown only for a limited region around the origin, while at 
infinity it approaches the unit value. Right: the charge and current distributions for this solution 
exhibits essentially the same structure as for the monopole rings, shown in FigJSJ 



Since there is no electric field for these solutions, their angular momentum fl4.66p is zero. 
However, in the limit where the Higgs potential vanishes, the Lagrangian (13.501) admits the 
global symmetry [27]| , 14( 



Aq -^ Aq cosh 7 + $ sinh 7, 



$ ^ $ cosh 7 + y4o sinh 7, 



A, 



A, 



(5.132) 



and this allows one to produce solutions with an electric field starting from purely magnetic 
solutions. As was noticed in 88|] , this transformation may also generate a non-zero angular 
momentum 

J ~ sinh 7, (5.133) 

provided that the original purely magnetic solution is not spherically symmetric and does not 
satisfy the first order Bogomol'nyi equations 27|. The latter condition excludes from consid- 
eration the (multi) monopole solutions of these equations. However, applying the symmetry 
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(15.1321) to the monopole-antimonopole solution, which does not fulfill the Bogomol'nyi equa- 
tions and whose Higgs field is long-range in the limit of vanishing Higgs potential, gives an 
electrically charged system, since both monopole and antimonopole then receive an electric 
charge of the same sign 85|| • The total electric charge calculated with the gauge-invariant 
definition (13.581) is Qci ~ sinh7, and applying the surface integral formula (14.661) gives the 
angular momentum 51 1 



J 



47rQe 
9 



(5.134) 



directed along the symmetry axis passing through the monopole and antimonopole. 
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FIG. 21: Schematic visualisation of the rotating monopole-antimonopole pair as superposition of 
a monopole-electric charge pair {+P,q) and an antimonopole-electric charge pair (— P, g), both of 
which making the same contribution to the angular momentum. 



It is interesting to understand the origin of this angular momentum. If it was an orbital 
angular momentum, then it would be orthogonal to the monopole-antimonopole symmetry 
axis, but it is parallel to it instead. Let us remember that a static system containing an 
electric charge q and a magnetic charge P has an angular momentum J = Pq independent of 
the distance between the charges and directed from g to P |138l]. This angular momentum can 
be obtained by integrating the Poynting vector as in Eg. (14.62 1 ) . We know that a rotating 
monopole-antimonopole pair contains two magnetic charges of the opposite sign and two 
electric charges of the same sign. It can therefore be visualized as a superposition of two 
pairs, (+P, g) and (— P, g), both of which making the same contribution to the angular 
momentum, as shown in Figj2Tl 

A similar construction of rotating solutions can be carried out using other known static 
solutions of the non-Abelian theory (I3.50p . as for example t he monopole rings described in 
Sec JIIIBl also for the generic Higgs field potential 1311], 105| . It turns out that the following 
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general relation holds 

J= '""^^-' [l + (-l)1 (5.135) 



where m,k are the winding numbers in the ansatz ( J3.56P 1311], 5l|, 105| . Therefore, 



using Eg. (13.59 1 ) ■ the angular momentum is zero in the presence of a nonvanishing magnetic 



charge, even though T° 7^ 105| . while the relation (I5.13ip is still valid for solutions with 
zero magnetic charge. 

Since they are related to the monopole-antimonopole solution known to be unstable 



161| . rotating monopole-antimonopole pairs are very probably unstable as well. However, 
the mechanism of their rotation is interesting and suggests, in particular, a possibility to 
have stationary, non-radiating spinning solitons without axial symmetry. 

Specifically, in theories without radiation, as in classical mechanics, steadily rotating 
bodies can be totally asymmetric. However, if they carry an electric charge distribution, 
say, then their electric dipole momentum will depend on time, hence producing radiation, 
unless they are axially symmetric. This suggests that non-radiating spinning field systems 
should be axially symmetric - the assumption usually made when studying spinning solitons. 



An interesting illustration of this has been found by Hen and Karliner 87|], who analyze 
situations where rapidly spinning systems can loose axial symmetry. A classical example 
of this phenomenon is provided by the rotating Jacobi ellipsoids in Newtonian gravity - a 
theory without radiation. Hen and Karliner study spinning baby skyrmions on a compact 
(sphere or disk) 2D space, in which case there is no radiation, and find that the system 
looses axial symmetry for large enough J. However, as was discussed above, in the limit of 
infinite M^, whe n ra diation can exist, spinning baby skyrmions remain axially symmetric for 



any value of J 133l | - since otherwise they would radiate. 

Now, the rotating monopole-antimonopole pairs suggest a possibility to have non- 
radiating, non-axially symmetric rotation even in theories with radiation. This can be 
achieved by simply taking three or more dyons not aligned along one direction. Of course, 
this should be realized by constructing a smooth, static, non-axially symmetric, finite energy 
solutions in a non-Abelian gauge field theory. If they exist, they may have a non-zero J due 
to the same mechanism as for the rotating monopole-antimonopole pairs. 
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VI. VORTONS 



We are now ready to explicitly construct vortons as localized, finite energy solutions of 
the elliptic field equations. By construction, they are stationary and non-radiating. Our 
results can be viewed as complementary to those obtained by Battye, Cooper an d Su tcliffe 



17( 1 (although in a completely different context) and by Lemperier and Shellard 118 1. 



A. The Witten model 



Vortons were originally suggeste d |46l| 



48 



45( 1 . 47( 1 in the context of Witten's model 



of superconducting cosmic strings 1771 ]. This model contains two Abelian vectors A 



(a) 



(a = 1, 2) interacting with two complex scalars (f) and a with the Lagrangian density 

^^ = ~\H FltJF^"^^" + D^^*D^^ + D^(j*D^(j - U. (6.136) 



4 



a=l,2 



.(a) 



|W 



|W 



Here F^J = df^Al — d^Ali are the Abelian field strengths, the gauge covariant derivatives 



i(i) 



of the scalars are -0^0 = {d^ — igiAli ')(j) and D^a = {d^ — ig2Aii ')a where gi and g2 are the 



(2)^ 



gauge coupling constants. The scalar field potential is 



u = -^x, 



vif + -^K\cr\\\a\^-2^i)+^\ 



?\a\\ 



(6.137) 



where A^, Ag- ,ri^,ri„ and 7 are positive constants. This theory is invariant under local 
U(l)xU(l) gauge transformations, so that there are two conserved Noether currents. 

The theory admits stationary, cylindrically symmetric solution s of the vortex t ype , sup - 



port ing a constant non-zero value of one of the Noether currents 1771], [13|, [86|, |9l|], [46|, 



4], [l32|. These are the superconducting cosmic strings. Vortons are supposed to be loops 
made of these strings. 

Within the full gauged model (I6.136P explicit vorton constructions have never been at- 
tempted - due to the complexity of the problem. However, the problem simplifies in the 
global limit of this model, for gi = g2 = 0. The gauge fields then decouple and the theory 
reduces to 

£ = df,(f)*df'(f) + d^a*d^a-U (6.138) 

so that the U(l) xU(l) internal symmetry becomes global. This global theory still keeps some 
essential features of the original local model. In particular, it still admits superconducting 
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vortex solutions, even though the corresponding current is now global and not local. One 
can therefor e stu dy global vortons made of these vortices, which was in fact the subject 
of Refs.|l7l|. (ll8l |. Below we shall construct the global vortons as stationary solutions of 
the field equations in the model (16.1381) . which has not been done before. The question of 
whether these solutions can be generalized within the full gauged model (16.1361) remains 
open. 

One can absorb t]^ in the definition of </>, a in (16.1381) to achieve t]^ = 1. Since the overall 
normalization of the potential could be changed by rescaling the spacetime coordinates, 
x'^ — > Ax'^, one can impose one more condition on the remaining four parameters A,^, X^, rif^, 7, 
although we do not use this option. 

A minimal value of the potential is achieved for |0| = t]^ = 1 and |cr| = 0, in which 
case U = 0. This minimum is global if 47^ > A^A^. The perturbative spectrum of field 
excitations around this vacuum consists of two massless Goldstone particles, corresponding 
to excitations of the phases of the fields, and of two Higgs bosons with the masses 



M^ = ^X^, M, = -y/7 - ^ X^7]l (6.139) 

The global U(l)xU(l) symmetry of the theory, (p -^ 0e'"^, a -^ cre*"^, leads to the 
conserved currents 

jf^) = 23fi(z0*9», jjt) = 23fJ(za*9^V), (6.140) 

with dfj_Jl\^ = d^J'fs = 0. The energy- momentum tensor is 

Tf,u = d^(j)*du(l) + du(l)*d^(j) + df,a*ducr + ducr*d^,a - g^^C , (6.141) 

where g^u is the spacetime metric. 

The Lagrangian field equations in the theory (I6.138P read 

^m5'^</>+^</> = 0, 9,9-a+^a = 0. (6.142) 

Expressing the complex scalar fields in terms of their amplitudes and phases as 

= /ie*'^\ a = /2e*'^^ (6.143) 

equations (I6.142p assume the form (a = 1,2) 

d.d^fa = {d,^ad>'^2)fa " ^ ^ ' d.Uld'^a) = 0. (6.144) 
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These equatio ns ad mit cylindrically symmetric vortex solutions discussed by several authors 



(see e.g. 46|], 118l |. 84l|). For these solutions one has 

= /i(p)e'"^ a = /2(p)e^(^^+'^*) , (6.145) 

where n G Z and p, u; G M, with ^ /i(p) -^ 1 and /2(0) ^ /2(p) ^ as <— p — »• cxd, 
respectively. The field (f) thus vanishes in the vortex core and approaches the finite vacuum 
value at infinity, while a develops a non-zero condensate value in the core and vanishes 
at infinity. The fields </>, a are sometimes referred to as vortex field and condensate field, 
respectively. The phase of changes by 2?™ after one revolution around the vortex, while 
the phase of a increases along the vortex. The z-dependence of the condensate field gives 
rise to a non-zero momentum along the vortex, P = J T^d?x ~ p. 

The qualitative vorton construction usually discussed in the literature suggests taking a 
piece of length L of the vortex and identifying its extremities to form a loop. The momentum 
along the vortex will then circulate along the loop, thus becoming an angular momentum, 
and the centrifugal force which arises is supposed to be able to compensate the tension of 
the loop, thereby producing an equilibrium configuration. The coordinate z along the vortex 
then becomes periodic and can be replaced by the azimuthal angle v?, so that p will have to 
assume discrete values, p = 271771/ L with m G Z. The central axis of the vortex where the 
field vanishes then becomes a circle of radius R = L/2'k. 

B. Vorton topology and boundary conditions 

The above considerations suggest describing the vortons by the ansatz 

= /i(p, z)e'^'^''^'\ a = f2{p, zy-^^^'-'. (6.146) 

The phase of a increases by 27r?Ti after one revolution around the z-axis, so that f2{p,z) 
should vanish at the axis for the field to be regular there. The phase ipi increases by 27rn 
after one revolution around the boundary C of the (p, z) half-plane, which is the 2;-axis plus 
the infinite semi-circle (see Fig j22l) . The regularity condition then implies that 0(p, z) must 
have n zeros somewhere inside C. Therefore any vorton configuration can be characterized 
by a pair of integers {n,m), which we shall call, respectively, vortex winding number and 
azimuthal winding number. 
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FIG. 22: Vorton topology. The phase of (p increases by Inn after one revolution around the 
boundary C of the half-plane, (p then vanishing n times somewhere inside C. In the simplest n = 1 
case the zero is located at a point {R, 0) corresponding to the vortex center. The vortex has a 
characteristic thickness a. 

It is instructive to compare these boundary conditions to those for the Faddeev-Hopf 
knots in Eg. (12.321) . They are the same - the phases of 0, a wind, respectively, along the 
orthogonal directions shown in Figl2l The vortons and knots have therefore similar topology. 
This suggests introducing the 'topological charge' analogues to the Hopf charge (12.121) 

N = nm. (6.147) 

Strictly speaking, this quantity will be topologically invariant, that is unchanged under 
arbitrary continues field deformations, only in the sigma model limit defined by the condition 
(16.1681) below. However, we shall call it topological charge in all cases and shall require the 
vortons to have A^ 7^ 0. 

Let us pass in ( 16.146P to spherical coordinates, 

= /i(r,^)e*'^^ ('■'''), a = f2ir,^)exp[i{mif + ujt)]. (6.148) 

Inserting this to Eqs. (l6.144p gives 



A/i 



Wif + ^(A^-i) + 7/n/i 



0, 



A/2 



771 



X 



r^ sin t9 



UJ 



+ -^Ul-vl) + ifl]h = 0, 



v(AV^i) = 0. 



(6.149) 
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At infinity the fields are required to approach the vacuum values, /i — ^ 1, /2 — ^ 0. One finds 
for large r 

^1 = ^:^ + . . . , /2 = - exp(-v/M2-cu2 r) + . . . , 



A^ C 

— (1 + 3 cos' ^) + ... + - 



/i = 1 + ^(1 + 3 cos' ^) + . . . + r exp(-M^ r) + . . . , (6.150) 



where A, B, C are integration constants and the dots denote the subleading terms. These 
asymptotic expansions show the presence of two massive Higgs modes and a long range 
massless Goldstone mode. The second Goldstone field is not excited within the ansatz 
06.1481) . Introducing the notation 

X = /icosV'i, F = /isinV'i, ^ = /2, (6.151) 

the equations assume the form 

'A 



A^ = ( V^ - ^' + ^(^' - ^l) + 7(^^ + Y^) ) Z. (6.152) 



j,2gjj^2^ 2 

Using (16.1501) . one has at large r 



X = l + 0(r-^), F = 0(r~'), Z ~ - exp(- ^M^ - cj2 r). (6.153) 

Introducing 

E2= f{\V(p\^ + \Va\^)£x, Eq= fud^x, M= f Z^£x, (6.154) 

the energy is 

E= T^d^x = u^Af + E2 + E0 , (6.155) 

while the Lagrangian 

L= I Cd^x = uj^Af-E2-Eo. (6.156) 

One of the two Noether charges is non- vanishing, 

Q = I d^xJ^^ = 2ioX, (6.157) 
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and this gives rise to the angular momentum 

J= T^d^x = mQ. (6.158) 

Rescahng the spatial coordinates as x ^ Ax, the Lagrangian L should be stationary, which 
implies the virial relation 

E2 = 3iuJ^^f-Eo). (6.159) 

The presence of the term (Tn^Z^/r^sin^-i?) in |Vo"p requires that for finite energy fields one 
should have Z = at the 2;-axis. Let us now remember that the phase ipi oi (p = X + iY = 
fic^^^ should increase by 2-7m after one revolution around the boundary of the (p, z) half- 
plane. Since one has X = l,y = Oatr^ oo, it follows that the phase is constant at 
infinity, so that it can only change along the 2-axis (see Fig j22!) . Therefore, 

ipi{p = Q,z = oo) — ipi{p = Q,z = — oo) = —27m, (6.160) 

from where it follows that the functions X and Y have (at least) 2n and 2n + 1 zeros at the 
2;-axis, respectively. Such a behavior is compatible with the assumption that X is symmetric 
and Y is antisymmetric under z — * —z. Assuming that Z is symmetric, one arrives at the 
following parity conditions, 

X{r,'d)=X{r,'K-d), r(r,^) = -F(r, vr - ^), Z(r,^) = Z(r, vr - ^), (6.161) 

which allow one to restrict in the analysis the range of d to the interval [0,7r/2]. Summa- 
rizing everything together, solutions of Eqs. (l6.152p should satisfy the following boundary 
conditions. 

At the symmetry axis, ^ = 0, functions X, Y should have, respectively, n and n — 1 zeros 
for < r < oo, and one should also have 

di)X = di)Y = Z = 0. (6.162) 

drX = Y = Z = 0. (6.163) 

X = 1, Y = Z = 0. (6.164) 

d^X = d^Z = Y = 0. (6.165) 



At the origin, r = 0, 

At infinity, r = oo. 

In the equatorial plane, 'd = 7r/2, 
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C. The sigma model limit 

We have succeeded for the first time to solve the equations in the special limit of the 
theory obtained by choosing the parameters in the potential f l6.137p as 

Aa = A^ = /3, ?7a = 1, l = -(3 + io, (6.166) 

such that the potential becomes 

^ = f(l0P + kP-l)' + 7o|0lVP. (6.167) 

Taking /5 — > oo enforces the constraint 

|0|2 + |(j|2 = i (6.168) 

and the theory ( 16.138^ becomes a sigma model, 

C = d^(j)*d''(j) + d^a*d>'a - 7o|0| Vr- (6-169) 

This model is simpler than the full theory (I6.138p . It has less degrees of freedom and no free 
parameters, since the value of 70 can be changed by rescaling the spacetime coordinates. 
The field mass M^- determined by (16.1391) reduces now to ^/to, while the mass M^ becomes 
infinite, which means that field degrees of freedom violating the constraint (I6.168P are ex- 
cluded from the dynamics. The theory therefore contains only one massive particle and two 
Goldstone bosons in the spectrum. It is worth noting that A^ = mn now becomes genuinely 
topological charge, invariant under arbitrary smooth field deformations. 

It is convenient to use the Lagrange multiplier method by adding to the Lagrangian 
(I6.169P a term /i(l — |0p — |crp) and varying with respect to (j),a and fi. Using the same 
ansatz (I6.148p . (l6.15ip as before gives 



AX = (70^2 + ^) X, 
AY = {-foZ^ + /i) Y, 

AZ = ^ . . -uj^ + 7o(X2 j^Y^) + ii]z. (6.170) 

and also 

X^ + F^ + Z^-l =0. (6.171) 
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Multiplying the three equations fl6.170p . respectively, by X, Y, Z and taking their sum and 
using the constraint (16.1711) one finds the expression for the Lagrange multiplier, 

2 

/i = XAX + YAY + ZAZ - 2-fo{X^ + Y^)Z^ - (- ^ - ij^)Z^ 

r^ sin 6 

- -A{X^ + Y^ + Z^), (6.172) 

where the second line has been added in order to cancel the second derivatives appearing in 
the first line. Inserting this to (I6.170p . the constraint (16.1711) will be imposed automatically 
on the solutions, so that it can be excluded from considerations from now on. 

The boundary conditions for Eqs. (l6.170l) are given by Eqs. (l6.162l) - (l6.165p . The regularity 
condition Z|£)=o,7r = imposes the constraint X"^ + F^ = 1 on the z-axis. In addition, since 
at the origin one has Y = 0, the constraint requires that X"^ = 1, and so one can choose 

X\^^o = -1 • (6.173) 

Since one has X\r=oD = +1, this will guarantee a non-zero value of n, which is the principal 
advantage of the sigma model theory (16.1691) . 

In fact, making sure that the phase of (p winds n ^ times around the contour C is 
not simple. This condition can be naturally implemented in toroidal coordinates (12.130 . 
but these coordinates are somewhat singular at infinity. Spherical or cylindrical coordinates 
are better suited for numerics, but the field topology is then determined by zeros of the 
X,Y amplitudes (see Eq. (16.160p ). whose number cannot be generically prescribed, since 
they are able to disappear during numerical iterations. Now, the condition (16.1730 enforces 
a non-trivial field topology at least in the n = 1 case. 

We can now solve the problem. Starting from a field configuration satisfying the boundary 
conditions (16. 1620 - 06. 165p . (16.1731) we numerically iterate it until the convergence is achieved. 
The resulting solutions are qualitatively very similar to the generic vortons described in detail 
below, and they also agree with the 'skyrmions' in the theory of Bose- Einstein condensation 
previously obtained in a completely different way by Battye, Cooper and Sutcliffe 17|. The 
latter issue will be discussed below in more detail. 

D. Explicit vorton solutions 

Having obtained vortons in the sigma model limit, we can use them as the starting profiles 
in order to iteratively descend from infinite to finite values of /5 in the potential (16.1670 , thus 
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relaxing the sigma model condition. After this, we can relax also the conditions (16.1661) . 
thereby recovering vorton solutions within the full original theory (16.1381) . 

In our numerics we use the program FIDISOL (written in Fortran), based on the iter- 
ative New ton- Raphson method. A detailed presentation of the FIDISOL code is given in 
152| , 1501 1 . In this approach the field equations are discretized on a (r, i}) grid with A^^ x N^ 
points, and the resulting system is iterated until convergence is achieved. The grid spacing 
in the r-direction is non-uniform, while the values of the grid points in the angular direction 
are given by dk = {k — 1)tx / {2{N^ — 1)). Instead of r, a new radial variable x is introduced 
which maps the semi-infinite region [0, oo) to the finite interval [0,1]. There are various 
possibilities for such a mapping, a flexible enough choice being x = r/{c + r), where c is 
a properly chosen constant. Typical grids for the n = 1 solutions have around 200 x 30 
points. The typical numerical error for the solutions is estimated to be of order 10"'^. In 
addition, the virial relation (I6.159P was also monitored. The deviation from unity of the 
ratio E2/3{uj'^N' — Eq) is less than 10~^ for most of the solutions we have considered. 




FIG. 23: Levels of constant amplitude \cj)\ (closed lines) and phase ipi (radial lines, red online) 
on the p (horizontal direction)-^ (vertical direction) plane for the vortex field (j) = \cj)\ exp(ii/;i) 
corresponding to the vorton solution with n = 1, m = 2, w = 0.85, A^ = 41.12, X^ = 40, rjfj = 1, 
7 = 22.3. 



As a result, for given values of the parameters A^, A^-, tj^j and 7 in the potential, we can 
specify the azimuthal winding number m and the frequency uj and obtain a vorton solution. 
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The value of n, the charge Q and energy E are computed from the numerical output. A 
complete analysis of the parameter space of solutions is a time consuming task that we 
did not aim at. Instead, we analyzed in detail a few particular classes of solutions, which 
hopefully reflects the essential properties of the general pattern. 

We mainly considered the case where the function X vanishes once at the positive z semi- 
axis, which corresponds to the vortex winding number n = 1. The amplitude of the scalar 
field (f) has in this case one zero located on a circle in the § = 7r/2 plane. The corresponding 
vorton topology can be illustrated by the diagram in Figj23| where the level lines of the 
complex vortex function = |0| expliipi) are shown. The levels of constant amplitude, |0|, 
are closed lines encircling the center of the vortex where (p vanishes. Emanating from this 
center there are lines of constant phase. The phase, ■i/^i, increases by 27r after one revolution 
around the center. 





FIG. 24: The amplitudes X{p, z),Y{p,z) for the typical vorton solution; here n = 1, m = 2, 
uj = 0.85, A^ = 41.12, \a = 40, r?^ = 1, 7 = 22.3. 



All vortons have toroidal distributions of the energy density and charge. It seems that 
they do not exist for arbitrary values of the parameters of the model but only for some regions 
in the parameter space. In Figs. [24 



I the 3D plots of X, F, Z and also of |0| = ^X^ + Y^ 
for a typical solution with m = 2 are presented. One can see that these functions exhibit 
a strong dependence both on p and z, with the zero of |0| located in the z = plane at 
i?~0.82. 



In Fig. [261 [27] the 2D profiles of X, Y,Z and |0| for a typical m = 1 solution for several 
values of ^ are shown; these pictures remain qualitatively the same for higher values of m. 
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FIG. 25: The amplitudes Z{p,z) and \(j){p,z)\ for the same solution as in Fig, 
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FIG. 26: Profiles of X{r,i!)),Y{r,'&) for fixed ?? for the n = m = 1 vorton solution with uj = 0.8, 
A^ = 33.26, A„ = 32.4, 77^ = 1, 7 = 18.4. 



As seen in these plots, the space can be partitioned into three regions. In the first region, 
located near the origin, the functions X, Y, Z present a strong variation, with Y approaching 
its extremum. The second region is located inside the vortex, where X stays very close to 
zero, while the field Z is almost constant and close to its maximal value. As seen in the 
insertion in the first plot, X crosses zero value with a non-zero first derivative, as it should, 
since the degree of zero of at the vortex center should be one. In the third region, outside 
the vortex ring, the fields approach the vacuum values. Both X and Z show strong variations 
in a transition domain between second and third regions. We notice also that the value of 
X at the origin, which was —1 in the sigma model limit, is now larger, although it stays 
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FIG. 27: Profiles of Z{r, -d), \(j){r, 'd)\ for fixed values of "i? for the same solution as in Fig, 



always negative and smaller that —0.5 for all solutions we have found. 





FIG. 28: The energy density Tq for the n = m = 1 (left) and n = 1, m = 3 (right) vortons with 
uj = 0.8, A^ = 39.25, A^ = 38.4, r/^ = 1, 7 = 21.39. 



For the m = 1 solutions the energy density is maximal at the origin. Emanating from this 
central maximum there is a toroidal shell of energy containing the ring of radius R where 
vanishes (see Fig l28|) . When m increases R grows, whereas the height of the maximum at 
the origin decreases (see Fig j28|) . so that the shape of the energy density resembles a hollow 
tube (see Fig j30|l . There is, however, an almost constant energy density inside the tube, 
whose value decreases with m. As seen in Figj29l the charge of the solutions, although also 
localized in a compact region, does not present the tube structure of the energy density. 

Nontrivial vorton solutions are likely to exist for all values of m. We have studied solutions 
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FIG. 29: The charge density for the same solutions as in Fig 
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FIG. 30: Left: the surface of constant (Tq = 1.05) energy density for the same n = 1, m = 3 vorton 
solution as in Figl281 It has the structure of a torus containing inside another toroidal surface. 
Right: the energy E{uj) and charge Q{u;) for the n = l,m = 2 vorton solutions. 

up to m = 5, but beyond this value the numerical accuracy of our procedure decreases 
drastically. The energy of the solutions increases with m, and for solutions with the same 
u one has E{m) < mE{m = 1). All solutions we have found have the ratio R/a between 
1 and 1.5, so that the radius of the vortex ring is almost the same as the thickness of the 
vortex. The construction of vortons in the thin ring limit, for R ^ a, remains a numerical 
challenge. 

Figj30] (right) shows the energy E and charge Q as functions of the frequency u for m = 2 
vortons. These functions are actually quite similar to those found in the Q-ball case. As 
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one can see, solutions exist for a limited range of frequencies, u^ < u < u^. Unfortunately, 
the numerical accuracy decreases considerably near the limits of this interval. However, it 
seems that both Q{uj) and E{uj) diverge in these limits, as in the Q-ball case. 

In analogy with Q-balls, the vortons become large as cu — i> uj±. However, it is not the 
vorton radius R that becomes large, since it appears to be not very sensitive to u and remains 
finite and nonzero when uj approaches its limiting values. It appears that the torus thickness 
increases for uj ^ uj^ so that its volume grows, thereby increasing the energy and charge. 
Both E and Q assume their minimal values aX uo = uOcnt- The dependence E{Q) appears to 
be double- valued, very similar to the one shown in FiglH] for Q-balls, again exhibiting two 
branches with a cusp. Solutions with oj < uj ait are less energetic than those for uj > uJcth 
and with the same Q. 

Using again the analogy with Q-balls, it seems likely that only solutions from the lower 
branch can be stable, probably those for large enough charge. Stable vortons certainly 
exist, perhaps not for all values of the parameters of the potentials, but at least for large 
enough values, when the theory approaches the sigma model limit. Vortons in the sigma 
model limit are stable, since, as will be discussed below, they can be obtained via 3D energy 



17j ] . This suggests that they should remain stable at least when they are close 



minimization 
to this limit. 

Although we did not study in detail the dependence of the solutions on the parameters 
of the potential, one should mention that for r]cr — 1 and uj ^ 0.85 (the values we mainly 
considered) we could find solutions only for A^ > 22, while the ratio \^/\a was always close 
to one. In other words, our vortons are not that far from the sigma model limit, which 
suggests that they could be stable. 

Vortons with n > 1 also exist, and we were able to construct them for n = 2. For these 
solutions there are two concentric rings in the d = 7r/2 plane where |0| vanishes. All n = 2 
solutions we have found have very large values of A^, X^, 7. 

Finally, an interesting class of solutions of Eqs.f l6.152l) is obtained by setting Y = 0, which 
means that the phase of is trivial and so n = 0. Such solutions have not been discussed 
in the literature. Although m 7^ in this case, one has A^ = nm = 0, and so according to 
our definition these solutions are not vortons, since they are 'topologically trivial'. However, 
they also have a ring structure, resembling somewhat spinning Q-balls. One can view them 
as '</)-dressed' Q-balls made of the complex a field, the non-renormalizable \a\^ term in the 
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Q-balls potential being replaced by the interaction with a real scalar field 0, as was first 
suggested in Ref. [75|. 



VII. RING SOLITONS IN NON-RELATIVISTIC SYSTEMS 



So far we have been considering relativistic ring solitons that could be relevant in models 
of high energy physics. At the same time, similar objects also exist in non-relativistic physics, 
and below we shall consider some of them which are most closely related to our previous 
discussion. 

A famous example of soliton-type toroidal systems is provided by the magnetically con- 
fined plasma in the TOKAMAK's. However, since considering this subject would lead us 
far beyond the scope of our pres ent discussion, we shall restrict ourselves to simply giving a 
reference to a recent review 281]. In addition, since the external field is necessary to confine 
the plasma in this case, this example does not quite fit in with our discussion, because we 
are considering closed, self-interacting systems. Another example of plasma confinement, 
although somewhat controversial, but which could perhaps be explained by some solitonic 
structures, is provided by the phenomenon of ball lightning (see 82| for a bibliography). It 
has in fact been conjectured that some analogs of Faddeev-Skyrme solitons in plasma might 
be responsible for their existence 67|], 63|], although we are unaware of any confirmations of 
this conjecture via constructing explicit solutions. 

Field theory models of condensed matter physics are more closely related to our discus- 
sion. As in the relativistic case, they can be either theories with local gauge invariance, as 
for example the Ginzburg-Landau theory of su per co ndu ctivity 771 ]. or models with global 
internal symmetries, as in the Gross-Pitaevskii 78|], 1351 ]. 83l ] theory of Bose-Einstein con- 
densation. It has been conjectured that some analogs of the Faddeev-Skyrme knot solitons 
could exist in the multicomponent Ginzburg-Landau models HI]. However, we are again 
unaware of any explicit solutions, and moreover, as was discussed above, their existence is 
not very plausible, at least in the purely magnetic case. On the other hand, global models 



of condensed matter physics do admit ring type solitons. 
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A. Vortons versus 'skyrmions' in Bose-Einstein condensates 

Surprisingly, it turns out that the described above vortons in the global limit of Witten's 
model can also be considered as solutions of the Gross-Pitaevskii (GP) equation, 



I- 




^5^K,fe|^,n^,. (7.174) 



dt \ 2 ' ' 2 

This equation, also called nonlinear Schrodinger equation, provides an effective description 
of the Bose-Einstein condensates. Here ^ a are the condensate order parameters, V^(x) is 
the external trapping potential and Hah describe interactions between different condensate 
components. In the case of a two-component condensate the indices a, h assume values 1, 2. 
Let 

(f) = 0(x), a = a(x)e*''* (7.175) 

be a solution of the equations (I6.142p of the global model of Witten (16.1381) in the case where 
the parameters in the potential (16.1371) are related as 

A^ = X.vl- (7.176) 



Then setting 



^1 = e-'- * 0(x), ^2 = e 'U + 2 ;* ^(x) (7.177) 



gives a solution of the GP equation (17.1741) in the case where 

V^(x) = 0, Kn = y, ^22 = y, Ki2 = K21 = 7. (7.178) 

It follows that any vorton solution satisfying the condition (17.1761) (this condition is typically 
not difficult to achieve numerically) provides a solutions to the Gross-Pitaevskii equation 
and hence can be interpreted in terms of condensed matter physics. 

In particular, the condition (17.1761) is achieved in the sigma model limit (I6.166p . (l6.168l) . 
in which case one obtains solutions of the GP equation with 

V{x) = Kn = K22 = 0, K12 = fi2i = 7o (7.179) 

by setting \l/i = 0(x), \l/2 = e~*^"*cr(x). Remarkably, solutions of the GP equation in this 
case have been independently studied by Battye, Cooper and Sutcliffe (BCS) |l7|, so that 
we can compare our results in Sec JVICl with theirs. They parametrize the fields as 

0(x) = cos f e'"'^ , a(x) = sin ^ e'"*^ , (7.180) 
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which exactly coincides with the parametrization fl2.32p of the Faddeev-Hopf field. Instead 
of directly solving the GP equation, BCS minimize the energy 

E2 + Eo= f{\V(P\^ + |Va|2 + 70101 Vn^^x (7.181) 

in a given A^ sector by keeping fixed the particle number 

^f= [\a\^£x. (7.182) 



They call the resulting configurations (f)j^{'x), cr^(x) 'skyrmions' - because Eg. (17. 1801) agrees 
with the skyrmion parametrization (I5.123p . At first view it is not completely obvious how 
these 'skyrmions' are related to our vortons 0a; (x), (Ta;(x) obtained by solving Eqs. (l6.170l) for 
a given u. However, these are actually the same solutions. In particular, the profiles of the 
vortons described in the previous sections, their energy density distributions, are similar to 
those for the 'skyrmions' given in Ref.|l7l|. 

The two different ways to obtain the solutions, either by minimizing the energy or via 
solving the equations, are actually equivalent^. This issue has in fact already been discussed 
above in the Q-ball context. Q-balls can be obtained either by solving Eq. (15.831) . which gives 
00; (x), or via minimizing the truncated energy functional Eq + E2 (15.891) with the particle 
number J\f (I5.88P fixed, which gives (j)j\j-{x.). Once one knows the value u = uj{JV), obtained 
either with the formula (15.901) or from the virial relation (I5.82p . the two results can be related 
to each other: (j)^{x) = (j)^(^j^){-x) = 0Ar(x). 

Similarly, the BCS skyrmions (f>j\j-{'x),aj\j-{'x) can be related to our vortons 0^; (x) , cr^; (x) 
by simply establishing the relation uj{N'). The knowledge of u is important to make sure 
that the no radiation condition 

cu^ <M^= 70 (7.183) 

is fulfilled, which is the case for our solutions, while their profiles agree with the features 



described by BCS in ITJ]. This provides an independent confirmation of our numerical 
results described in Sec JVI C\ Since they can be obtained by the 3D energy minimization, 
these solutions are stable. 



Apart fro m th e work o 



BCS ther e hav e been other studies of 'skyrmions' in Bose-Einstein 



condensates [149|, [147|, |l44| . (991], |l24[ |. Such solutions have been constructed by numeri- 



cally resolving the GP equation for more general choices of V^(x) and Kab and also for more 



^ We thank Richard Battyc for explaining this point to us. 
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than two condensate components. Possi biliti es for an experimental creation and observation 
of such objects have also been discussed |l46| . |l45| . In all studies the solu tions are typically 
presented for the winding numbers n = 1, with the exception of Ref. 144l |. where results for 
n = 2 are reported, although without giving many details. In our language this corresponds 
to the double vortons. 



B. Spinning rings in non-linear optics — Q-balls as light bullets 



Yet another application of the considered above ring solitons arises in the theory of light 
propagating in media whose polarization vector V depends non-linearly on the electric field 
S. The sourceless Maxwell equations become in this case non-linear and in a number of 
important cases they can be reduced to a non-linear Schrodinger (NLS) equation similar to 
the Gross-Pitaevskii equation (17.1741) . For a Kerr medium, when the non-linear part of V 
is cubic in S, the NLS equation has exactly the same structure as Eq.f l7.174|) . but for more 
general media it can have different non-linearities. Soliton solutions of this NLS equation 
describe non-linear light pulses, sometimes called light bullets, which are very interesting 
from the purely theoretical viewpoint and which can actually be observed, as for example in 
the optical fibers. Unfortunately, discussing these issues in more detail would lead us away 
from our subjects, and so we simply refer to a monograph [o] on optical solitons. 

Most of the known solitons of the NLS equation describe plane waves or cylindrical beams 
of light 6| . However, solutions describing solitons l ocali zed in all three dimensions are also 
known, they have been studied by Mihalache et al 125|. These solutions describe spinning 
rings, which is very interesting in the context of our discussion. 

Mihalache et al consider light pulses travelling in the z direction with a group velocity 
y in a medium with cubic and quintic non-linearities. After a suitable rescaling of the 
coordinates and fields, the envelope of the pulse, \1/, satisfies the equation 125 1 



dz 



d^ 



d^ 



d^ 



dx"^ dy'^ dC/ 



+ 1^1 



|*P * 



(7.184) 



where C, = z — Vt. Solutions of this equation, numerically constructed in Ref. 125[ |. describe 
'spinning light bullets'. Interestingly, it turns out that these solitons correspond to the 
spinning Q-balls already described above. 

Specifically, let us rename the variables a.s z ^>- 1, ( ^ z and set \1/ = exp{z(u;^ — M'^)t}(p. 
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Eq. fl7.184p then assumes the form 

i0 = (-A + 101^ - |0p + M^ - uj^)(j) . (7.185) 

Let us now consider Eq. fl5.73p describing Q-balls with the potential (15.761) . We notice that 
rescahng the spacetime coordinates and the field as x^ —>■ Aix'^, $ — »■ A2$ is equivalent to 
changing the values of the parameters A, a in the potential. Therefore, as long as they do 
not vanish, one can choose without any loss of generality A = 1/3, a = 3/2. Setting then 
$ = e*'^*0, Eq. (15.731) assumes the form 

- (^ - 2iu^ = (-A + 101^ - |0|2 + M^ - uj^)(f) , (7.186) 

with M^ = Xb. We see that if = then Eq. fl7.185p and Eq. fl7.186p become identical. They 
have therefore the same solutions. We thus conclude that the Q-balls solutions discussed 
abo ve in Sec JVAl can describe light pulses in non-linear media. The solutions found in 
Ref. 125| correspond to the simplest m~^, that is non-twisted, even parity Q-balls. Let us 
remind that for a given m these solutions can be parametrized by the value of u varying 
within the range u;^(?7i) < u'^ < u"^ = M^. 

Mihalache et al also study the dynamical stability of the solutions by analysing their 



small perturbations 125l |. Perturbing in Eq. (17. 1851) the field as (p ^ (p + 6(f) and linearizing 
with respect to 6(j) gives 

i6(j) = V[(f)]S(j) , (7.187) 

where 'D[4>] is a linear differential operator. If the spectrum of this operator is known, 

T)mx = \^x, (7.188) 

then Eq. (17. 1871) is solved by setting S(f> = exp{iAt}^A- It follows that if there are complex 
eigenvalues A with a negative imaginary part, then the perturbations grow for t ^ oo and 
so the background is unstable. Mihalache et al study the spectral problem (I7.188P and find 
complex eigenvalues for all spinning solitons with m = 2. They also integrate Eq. (l7.184p to 
trace the full non-linear perturbation dynamics, and they observe that the spinning toroidal 



soliton splits into three individual non-spinning spheroidal pieces [125| . However, they find 
that the ?n = 0, 1 solitons can be stable, both at the linear and non-linear levels, if u is 
less than a certain value u;*(rra) < M but still within the domain of the solution existence, 
u;_(m) < a; < uj^{m). According to their results, for m = one has u^ = cUcrit, which is the 
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value where E{uj),Q{uj) pass through the minimum (see Fig JT^ . For m = 1 the stabihty 
region shrinks, since uj^ < tUcrit, so that not all parts of the less energetic solution branch 
(see FiglHl) correspond to stable solutions. 

These results are very interesting, since they allow us to draw certain conclusions about 
stability of the relativistic Q-balls. Linearizing Eq. (I7.186P gives 

-50-2zcj(50 = P[0]50, (7.189) 



whose solution is 5(f) = exp{i'~ft}^x with 7 = —uj ± VA + cJ^. AH eigenmodes with complex 
A therefore correspond to complex 7 and so to unstable modes. In addition, there could 
be unstable modes with real A < — tu^, since 7 would then be complex. As a result, the 
relativistic Q-balls have the same unstable modes as their non-relativistic optical counter- 
parts (although with different instability growth rates) and perhaps also some additional 
instabilities. 

It follows that the 2"*" spinning Q-balls are unstable, while the l"*" ones, if they belong 
to the lower solution branch and their charge is large enough, should be stable (unless they 
have negative modes with A < — u;^). The latter conclusion agrees with the results shown in 
FiglH 

As was said above, Q-balls corresponding to the lower solution branch in FiglH] cannot 
decay into free particles. As we now see, this property does not yet guarantee their stability, 
since there can be other decay channels. Their non-relativistic optical counterparts can decay 
into several solitonic constituents. In order to see if the same could be true for the relativistic 
Q-balls, one should integrate the full time evolutions equation Eq. (17. 1861) . The existence of 



negative modes for sp innin g Q-balls is also suggested by t 



he results of Ref.[l0(,ll6l. 



Q: 



126|,[127| to describe the case of 



The results of Ref . |l25l | have been generalized in Refs. p 
light pulses with two independent polarizations and also for more general non-linearities of 
the medium. 

C. Moving vortex rings 

For the sake of completeness, we would also like to consider another well known type of 
ring solitons: moving vortex rings in continuous media. Although they are stabilized by the 
interaction with the medium and not by intrinsic forces, their structure can be quite similar 
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to that for the other sohtons considered above. 

A vortex ring moving in a medium encounters the 'wind' that produces a circulation 
around the vortex core, giving rise to the Magnus force - the same force that acts on a spin- 
ning ping-pong ball in the direction orthogonal to its velocity. This force is directed outwards 
orthogonally to the ring velocity, stabilizing the ring against shrinking (see Fig l3Tl) . As a 
result, the ring travels with a constant speed (if there is no dissipation), dynamically keeping 
a constant radius, so that in the comoving reference frame its configuration is stationary. 



moving vortex ring 





vorton 



FIG. 31: Moving vortex rings are stabilized by an analog Fm of the Magnus force produced by the 
phase rotation around the vortex core. Vortons are stabilized by the centrifugal force produced by 
the phase rotation along the vortex. 



The moving vortex rings are best known in fluid dynamics, as for example smoke rings. 
The vortices themselves correspond in this case to their traditional definition - these are line 
type objects containing a non-zero rotation of the fluid in their core, as for example swirls 
in water or tornadoes. Vortices and vortex rings can be created in the wake of rigid bodies 
travelling through the fluid. Their hydrodynamical theory was created in the 19-th century 
by Helmholtz and Kelvin and has been developing ever since, finding numerous practical 
applications, as for example in the aircraft engineering. We do not intend to discuss this 
vast subject here and simply refer to the standard monograph on vortices 
give the classical formulas of Kelvin relating the vortex ring velocity V and its energy E 
and momentum P to the ring radius, R, and the vortex core radius, a, 

r 



148|. We also 



2 ' V a 4 



P = Ti'qTR-' 



AttR \ a 4 



(7.190) 



which are valid in the thin ring limit, R ^ a. Here rj is the fluid mass density and F is 
the circulation of the fluid velocity around the vortex. These formulas show that large rings 
have large energy and momentum but move slowly. Very large and slowly moving rings, 
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which could perhaps be described by some version of hydrodynamics of magnetized plasma, 
are observed in Sun's atmosphere. 

1. Moving vortex rings in the super fluid helium 



Moving vorte x rin gs ex ist also in the superfluid helium, where they are created by moving 
impurities 140l | (see 17Cll | for a survey of defects in helium). It turns out that large rings 
in helium can be well described by the classical hydrodynamics. However, in the generic 
case their adequate description should be quantum. The superfluid helium can be approx- 
imately modeled by a weakly interacting Bose-Einstein condensate (although in reality the 
interactions in helium are not weak, see e.g. 23|), in which case vortices can be described 
by solutions of the Gross-Pitaevskii equation (17.1741) . In the simplest case one can consider 
the one-component version of this equation, which, upon setting V^(x) = 0, kh = n and 
^1 = e~*2* ^ reduces to 

(7.191) 



2z— =(-A + k(|^|2-1))*. 
Rescaling coordinates and time one can set ft = 1, if it is positive. Choosing ^ = e*"'''/(p) 



gives 



P 



n 



/ = (/'-!)/■ 



(7.192) 
1 describes a global 



whose solution with the boundary conditions /(O) = and f{oo) 
vortex with n quanta of circulation 78 1. 

Solutions of the GP equation (I7.19ip describing moving vortex rings and their various 
properties have been comprehensively analyzed by P.H. Roberts and collaborators in a series 
of articles written over the last 25 years, of which we only mention Refs. Sjj, 97|], 2J|, 22| . 



2l|. The fundamental vortex ring solutions have been numerically obtained by Jones and 



Roberts 



97( 1 by setting 



^ = ^(p,z-1/t), 
such that \l/(p, z) satisfies the equation 

- t2Vd,'i' = (-A + (|^|2 - 1)) ^ • 



(7.193) 



(7.194) 



They imposed the reflection condition \E'(p, z) = \E'(p, —z:)* and also the asymptotic condi- 
tions for large r, 

3fJ(^) = l + 0(r-3), 3(^)=0(r-2). (7.195) 
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Numerically integrating Eg. (17. 1941) . they obtained a family of solutions labeled hj V & 
(0, Knax), where V^max = l/"\/2 is the speed of sound in the dimensionless units chosen. For 
small V the solutions have ring profiles: the function \1' vanishes at a point [R, 0) in the 
(p, z) plane while its phase increases by 27r after one revolution around this point, so that \E' 
has the 'vorton topology' (see Fig l22l) . When V is small, the radius of the ring, R, as well 
as its energy and momentum, 

P^ = Q /(^* - 1)9^^ d^x, (7.196) 



are large, and in the leading order they agree with the classical Kelvin formulas (17. 1901) 81|. 
As V increases, E, P^ and R all decrease, and the following relation holds, 

V-gj,. (7.197) 

The ring radius R vanishes for V = Vq ^ 0.62 and for Vq < V^ < Vm^x solutions change their 



character and describe acoustic excitations in helium, called in 97| rarefaction pulses, in 
which case \E' no longer has zeroes . The energy and momentum diverge when V approaches 
0, Knax, and they both assume their minimal values for some V = V^rit ~ 0.657. The 
functions E{V) and PsiV) exhibit therefore qualitatively the same behaviour as E{uj) and 
Q{uj) for the Q-balls (see Fig JT^ . The function E{P^) shows the two-branch structure with 
a cusp, similarly to E{Q) in FigJT4l The moving rings belong to the lower branch, which 
suggests that they could be stable. Their perturbative stability has been demonstrated by 



Berloff and Roberts 2J], who also studied the r ing creation by solving the GP equation for 



a superfiuid flow around a moving solid sphere 22| and showed that rings are created in its 
wake, provided that its velocity exceeds some critical value. 

The moving ring solutions have been generalized for a non-zero trapping potential ^(x) 



in the GP equation (see |110| for a review) and for the two-component condensate case 



[20| . An interesting solution of the GP equation describing a non-linear superposition of a 
straight vortex and a vortex ring encircling it and moving along it, somewhat similar to the 



'hoop' solution shown in FigJT9l has been presented by Berloff in 21|. 

One may also wonder whether the vortex ring solutions could be generalized within a 
relativistic fleld theory (as is the case for the vortex solution of Eq. (17. 1921) ). In fact, the 
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answer to this question is affirmative, since taking a solution \l/(p, z) of Eg. (17. 1941) and setting 

= e'^("^*)^(p,2) (7.198) 

and also a = solves Eqs. (l6.142l) of Witten's model for A^ = 2. However, since \E^ — i> 1 at 
infinity, it follows that (p -^ Q^y{z±t) -^^ ^j-^jg \\Ya\i^ which corresponds to a massless Goldstone 
radiation. The relativistic analogs of the moving vortex rings are therefore radiative solutions 
and their total energy is infinite. 

2. Moving magnetic rings 



An in teres ting example of moving ring solitons, found by Cooper 4l| and studied also by 



Sutcliffe 159l | , arises in ferromagnetic systems whose dynamics is described in the continuum 
limit by the Landau-Lifshitz equation 

dn°- 



eabcn''An''. (7.199) 



dt 

Here n'^ is a unit three- vector, n""n"- = 1, giving the local orientation of the magnetization. 
The temporal evolution defined by the Landau-Lifshitz equation preserves the value of the 
energy 

E=jf{dkn''f£x. (7.200) 

In addition, it does not change the momentum. 

Pi = ^ f e^jkixjBk)d\ (7.201) 

where Bi = \eijkTjk and Tjk is defined as in the Faddeev-Skyrme model by Eq. (12.21) . The 
number of spin flips. 



^--2 



- f{l-n^)d^x, (7.202) 

is also preserved by the evolution; for this number to be finite one should have n'^ = 1 at 
infinity. Since the theory contains a unit vector field, it can be characterized by the topo- 
logical index of Hopf defined by the same equation (12.91) as in the Faddeev-Skyrme model. 
An axially symmetric field configuration with a non-zero Hopf charge can be parametrized 
exactly as in Eq. (l2.10l) (with n = 1), 

n^ + in^ = sin Q{p, z)exp{i[mip — ip{p,z)]}, n^ = cosQ{p, z). (7.203) 
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The energy density then has a ring structure and the Hopf charge is N = m. Using this to 
calculate the momentum fl7.20ip shows that its ^-component does not vanish, 

Pz = ^J p{d,ed,ij - d,ed,ij)d^^. (7.204) 

This means that if one uses (17.2031) as the initial data for the Landau-Lifshitz equation, 
;hen this data set will have a non-zero initial momentum. To verify this assertion, Sutcliffe 



159| directly integrates the Landau-Lifshitz equation to reconstruct the temporal dynamics 
of the ring. For largely arbitrary functions 0(p, z), ip{p, z) in (17.2031) he then sees the whole 
ring configuration drift with an almost constant velocity along the 2;-axis. He also observes 
that the phase of n^ + in^ increases as a linear function of the distance travelled by the ring, 
so that there are internal rotations in the ring as it moves. 

These empirical features can be understood by making the following ansatz for the dy- 



namical fields 411], 



n^+in^ = smQ{p,z-Vt)exp{i[mip + ujt-ip{p,z-Vt)]}, n^ = cosQ{p, z-Vt), (7.205) 

with constant V,uj. Inserting this ansatz to the Landau-Lifshitz equation (I7.199p . the t,ip 
variables separate and the problem reduces to the two dimensional elliptic problem for 

Q{p,z),ij{p,z), 



Ae = {uj + Vd,il)) sin 9 + 



sin 9 cos 9, 



P 
sin9AV^ = -1/(9,9 -2cosQ[dpQdpi) + d,Qd,%l)]. (7.206) 

This problem is in fact very much similar to that in the vorton case. This can be seen 
by applying the same change of variables which was used to rewrite the Faddeev-Skyrme 
model (12.11) in the CP^ form (12.291) . n" = (^'^t°-(^ , with $t$ = 1. This parametrization of 
the system allows one to write down a simple Lagrangian for the Landau-Lifshitz equation 

dnnnD, 

£ = i^^dt^ - -{dk-rff + /i($^<l> - 1), (7.207) 

where p is the Lagrange multiplier. One can parametrize $ exactly in the same way as in 

the vorton case, 

X + iY 

y piiot+irrup 
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Inserting this to the Lagrangian (17.2071) and varying with respect to X, Y, Z, fi gives equations 
similar to Eqs. fl6.170l) for vortons in the sigma model limit. These equations can be solved 
with the same boundary conditions fl6.162p - fl6.165l) for X, Y, Z as for vortons. On the other 
hand, setting X + iY = cos(6/2)e*'^ and Z = sin(6/2), as in Eq. (17. 1801) . gives again 
Eqs.(E2Q6l). 

Fixing m one can solve the equations for suitable input values of u, V. Having obtained 
the solutions, one can calculate E{u,V), A/'(u;, \^) and P^lujV). Alternatively, instead of 
solving the equations, one can minimize the energy E wit h fix ed A/", P3. Having found the 



minimal value E{Af, P3), one can reconstruct u, V as |4l|. [l59l | 



--W' ^-I- '--' 



Cooper 41j and Sutcliffe 1591 ] choose the latter approach: they minimize the energy with 
fixed A/", P3 for several values oi m = N. Their numerics converge to non-trivial configu- 
rations Q{p, z),iIj{p, z), inserting which to (17.2051) gives dynamical rings travelling with a 
constant velocity. Fixing all other parameters, the ring velocity decreases with growing m. 
In fact, these rings also have an angular momentum, J = mj\f, but apparently this 
feature is not so essential for their stabilization as it is for vortons, since solutions with 
m = N = J = also exist and do not exhibit any particular features as compared to the 



solutions with m ^ |l59l |. Similarly to the other moving ring solitons discussed above, 
magnetic rings are stabilized by a Magnus-type force produced by the phase circulation 
around the vortex core and directed outwards, orthogonally to the velocity (see Fig JSTIl . 
The magnetic rings are stable not only against shrinking but also with respect to all other 



deformations - since they can be obtained via the 3D energy minimization jl59|. 



VIII. CONCLUDING REMARKS 

We have reviewed the known field theory solutions describing stationary vortex loops. 
From the physical point of view, loops of magnetic flux would perhaps be the most in- 
teresting. However, as we have seen, almost all explicit solutions describe vortex loops in 
global field theories, which are effective theories with a relatively limited range of applica- 
bility. Perspectives for generalizing these solutions in the context of gauge field theory do 
not look very promising at the moment, at least as far as knot solitons are concerned. The 
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best known static knot solitons, Faddeev-Skyrme knots, do not seem to admit immediate 
gauge field theory generalizations, unless one makes additional physical assumptions to fix 
the charges in the Protogenov-Verbus formula. Such assumptions can be made, as shows 
the example of Schmid and Shaposhnikov. However this requires a rather exotic physical 
environment. 

The non-Abelian monopole and sphaleron rings are very interesting theoretically, how- 
ever, they are presumably unstable. It is therefore unclear whether they can find important 
physical applications. Summarizing, the existence of knots - magnetic vortex loops sta- 
bilized against shrinking by internal stresses - does not seem to be very plausible in the 
context of physically interesting relativistic gauge field theories. 

Perhaps gauged vortons have more chances to exist. As we have seen, their global coun- 
terparts do exist as stationary, non-radiating field theory objects. Already for these global 
vortons an additional analysis to study the structure of the parameter space is needed. A 
natural problem to address would be their analysis in the thin ring limit, for large values of 
the azimuthal winding number m, when the ring radius is much larger than the vortex core 
thickness. One can expect that the field theory solutions in this limit should agree with the 
effective macroscopic description. This limit, however, is difficult to explore numerically, 
since the vorton fields are then almost everywhere constant, except for a narrow ring region 
containing all the field gradients. It is difficult to properly adjust the grid in this case, so 
that some other numerical methods are necessary, as perhaps the multi-domain spectral 



method [80|]. Is seems also that some other methods are necessary to study solutions with 
n> 1. 

Yet another interesting related problem would be to find vortons with three winding 
numbers. We know that the phase of </> winds around the vortex core, while that of a winds 
along the ring, the corresponding winding numbers being n and m. However, nothing forbids 
the phase of a to wind both around the core and along the ring, as does the phase of the 
twisted Q-balls, in which case one would need a third integer, k, counting the windings of a 
around the vortex core. Constructing vortons with k ^ remains a challenge. 

A very important open problem related to all vortons is their dynamical stability analysis. 
This problem can be handled either by studying the spectrum of small field fiuctuations 
around the vorton background, or by looking for vorton negative modes within the energy 
minimization method, or by reconstructing the full temporal evolution of the perturbed 
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vorton configuration. Stable vortons definitely exist, at least in the sigma model limit, but 
in the generic case their stability should be studied. 

Generalizing global vortons within gauge field theory does not seem to be impossible. 
In theories where all the internal symmetries are gauged, as in Witten's model (16.1361) . all 
spinning phases of the fields can be gauged away, so that if the fields are stationary and 
axisymmetric, these symmetries are manifest. In this case the surface integral formula fl4.66p 
for the angular momentum applies, providing higher chances for J to vanish due to the fast 
asymptotic falloff of the fields. Although the example of spinning gauged Q-balls shows that 
spinning is possible in local theories, it seems that vortons have better chances to exist in 
models with both local and global internal symmetries. Such models contain gauge fields, 
such that one can have magnetic vortices, but the gauge freedom is not enough to gauge 
away all the phases of the scalars. The symmetries of the fields are therefore not manifest 
and the surface integral formula does not apply, in which case the angular momentum should 
normally be non-zero. 

For example, it would be natural to analyze the existence of stationary vortons in a 
half gauged version of Witten's modes obtained by coupling the two scalars in Eq.f l6.138|) 
to a U(l) gauge field. Such solutions could then perhaps be generalized within the fully 
gauged U(l)xU(l) Witten's modes to confirm the existence of stationary loops made of 
superconducting cosmic strings - the idea first put forward more than 20 years ago. 

However, a real challenge would be to construct stationary superconducting loops in 
Standard Model. There are indications that this might be possible. The electroweak sector 
of St andard Model contains stationary current carrying vortices - superconducting strings 



166| . Their stability has been studied so far only in the semilocal limit, where the SU(2) 



gauge field decouples and the current becomes global [73l|. |72[| . and it has been found that 
sufficiently short pieces of strings are perturbatively stable |76|. This suggests that small 
loops made of these strings may also be stable. If such loops could be constructed for the 
physical value of the weak mixing angle, this would give stable solitons in Standard Model. 
For the sake of completeness, having in mind possible applications of vortex loops in 
astrophysics and cosmology, it is also interesting to mention the effects of gravity on these 
solutions. On general grounds, one expects all known ring solitons to admit self-gravitating 
generalizations with essentially the same properties as in the zero gravity limit, at least for 
small enough values of the gravitational coupling constant. However, the soliton structure 
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can change considerably in the strong gravity regime |l67j ]. The sohtons can then become 
gravitationally closed and may also contain a small black hole in their center. The spinning of 
the latter can endow the whole configuration with an a ngular momentum, even if the soliton 



itself cannot spin - as for example the monopole |104| . However, if the soliton can spin on 
its own, as for example a vorton, it would be interesting to put a spinning black hole in its 
center. Since there are no asymptotically fiat toroidal black holes in four dimensions J96| . 
the resulting configuration is expected to be a black hole with spherical horizon topology, 
surrounded by a vortex ring. 
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